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, The physics during the inflationary stage of the universe is of quantum na- 

ture involving extremely high energy densities. Moreover, it is out of equilib- 
rium on a fastly expanding dynamical geometry. We present in these lectures 
non-perturbative out of equilibrium field theoretical methods in cosmological 
universes. We then study the non- linear dynamics of quantum fields in matter 
t^- [ and radiation dominated FRW and de Sitter universes. For a variety of initial 

conditions, we compute the evolution of the quantum inflaton field, its large 
<^ ' quantum fluctuations and the equation of state. We investigate the explosive 

■ particle production due to spinodal instabilities and parametric amplification 

in FRW and de Sitter universes with and without symmetry breaking. We 
show how the particle production is sensitive to the expansion of the universe. 
In the large N limit for symmetry breaking scenarios, we determine generic 
late time fields behavior for any flat FRW and de Sitter cosmology. We show 
that the amplitude of the quantum fluctuations fall off in FRW with the 

03 

square of the scale factor while the order parameter approaches a minimum 
• i— i ! of the potential in the same manner. We present a complete and numeri- 

cally accessible renormalization scheme for the equation of motion and the 
energy momentum tensor in flat cosmologies. Furthermore, we consider an 
O(N) inflaton model coupled self-consistently to gravity in the semiclassical 
approximation, where the field is subject to 'new inflation' type initial condi- 
tions. We study the dynamics self-consistently and non-perturbatively 
with non-equilibrium field theory methods in the large N limit. We find that 
spinodal instabilities drive the growth of non-perturbatively large quantum 
fluctuations which shut off the inflationary growth of the scale factor. We 
find that a very specific combination of these large quantum fluctuations plus 
the inflaton zero mode assemble into a new effective field. This new field 
behaves classically and it is the object which actually rolls down. The metric 
perturbations during inflation are computed using this effective field and the 
Bardeen variable for superhorizon modes during inflation. We compute the 
amplitude and index for the spectrum of scalar density and tensor perturba- 
tions and argue that in all models of this type the spinodal instabilities are 
responsible for a 'red' spectrum of primordial scalar density perturbations. 
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I. INTRODUCTION AND MOTIVATION 



Inflationary cosmology has come of age. From its beginnings as a solution to the hori- 
zon, flatness, entropy and monopole problems [[[]], it has grown into the main contender for 
the explanation of the source of primordial fluctuations giving rise to large scale structure 
26| . There is evidence from the measurements of temperature anisotropies in the cosmic 



microwave background radiation (CMBR) that the scale invariant power spectrum predicted 
by generic inflationary models is consistent with observations |2|-|J and we can expect fur- 
ther and more exacting tests of the predictions of inflation when the MAP and PLANCK 
missions are flown. In particular, if the fluctuations that are responsible for the temper- 
ature anisotropies of the CMB truly originate from quantum fluctuations during inflation, 
determinations of the spectrum of scalar and tensor perturbations will constrain inflationary 
models based on particle physics scenarios and probably will validate or rule out specific 
proposals P,|2"7fl. 

The tasks for inflationary universe researchers are then two-fold. First, models of in- 
flation must be constructed on the basis of a realistic particle physics model. This is in 
contrast to the current situation where most, if not all acceptable inflationary models are 
ad-hoc in nature, with fields and potentials put in for the sole purpose of generating an 
inflationary epoch. Second, and equally important, the quantum dynamics of inflation must 
be understood. This is extremely important, especially in light of the fact that it is exactly 
this quantum behavior that is supposed to give rise to the primordial metric perturbations 
which presumably have imprinted themselves in the CMBR. This latter problem is the focus 
of this review. 

The inflaton must be treated as a non- equilibrium quantum field . The simplest way to 
see this comes from the requirement of having small enough metric perturbation amplitudes 
which in turn requires that the quartic self coupling A of the inflaton be extremely small, 
typically of order ~ 10~ 12 . Such a small coupling cannot establish local thermodynamic 
equilibrium (LTE) for all field modes; typically the long wavelength modes will respond 
too slowly to be able to enter LTE. In fact, the superhorizon sized modes will be out of 
the region of causal contact and cannot thermalize. We see then that if we want to gain 
a deeper understanding of inflation, non-equilibrium tools must be developed. Such tools 
exist and have now been developed to the point that they can give quantitative answers to 



these questions in cosmology - [fig] , |T^j20|j32jj33|1 . These methods permit us to follow the 
dynamics of quantum fields in situations where the energy density is non-perturbatively 
large (~ 1/A). That is, they allow the computation of the time evolution of non-stationary 
states and of non-thermal density matrices. 

Our programme on non-equilibrium dynamics of quantum field theory, started in 1992 
0, is naturally poised to provide a framework to study these problems. The larger goal 
of the program is to study the dynamics of non-equilibrium processes from a fundamental 
field-theoretical description, by solving the dynamical equations of motion of the underlying 
four dimensional quantum field theory for physically relevant problems: the early universe 
dynamics, high energy particle collisions, phase transitions out of equilibrium, symmetry 
breaking and dissipative processes. 

The focus of our work is to describe the quantum field dynamics when the energy density 
is high. That is, a large number of particles per volume m~ 3 , where m is the typical mass 
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scale in the theory. Usual S-matrix calculations apply in the opposite limit of low energy 
density and since they only provide information on in — > out matrix elements, are unsuitable 
for calculations of expectation values. 

In high energy density situations such as in the early universe, the particle propagator (or 
Green function) depends on the particle distribution in momenta in a nontrivial way. This 
makes the quantum dynamics intrinsically nonlinear and calls to the use of self-consistent 
non-perturbative approaches as the large N limit, Hartree and self-consistent one- loop ap- 
proximations. 

There are basically three different levels to study the early universe dynamics: 

1. To work out the nonlinear dynamics of quantum fields in Minkowski spacetime. By 
non-linear dynamics we understand to solve the quantum equations of motion includ- 
ing the quantum back-reaction quantitatively @ - U, [|l^,^4|,[^,[33|] . This level is in 



fact appropriate to describe high energy particle collisions . 

To work out the nonlinear dynamics of quantum fields in fixed cosmological back- 
grounds [T^,[T3|. New phenomena arise then compared with 1. showing that a 
Minkowski analysis is not quantitatively precise for expanding universes. 

A self-consistent treatment of the quantum fields and the cosmological background 



ll.Q. That is, the metric is obtained dynamically from the quantum fields (matter 



source) propagating in the that metric. 

We shall successively present the three levels of study. The first stage was reviewed 
in the 1996 Chalonge School ||. The second level is the subject of sees. VI and VII. 
We study the parametric and spinodal resonances both in FRW and de Sitter backgrounds 
wide range of initial conditions both in FRW and de Sitter backgrounds [TT||T^] • [Parametric 



resonance appears in chaotic inflationary scenarios for unbroken symmetry whereas spinodal 
unstabilities show up in new inflation scenarios with broken symmetry]. Both types of 
unstabilities shut-off through the non-linear quantum evolution as described in sees. VI 



and VII [|lT| , |13| both analytically and numerically. We follow the equation of state of the 
quantum matter during the evolution and analyze its properties. 

The third stage of our approach is to apply non-equilibrium quantum field theory tech- 
niques to the situation of a scalar field coupled to semiclassical gravity, where the source of 
the gravitational field is the expectation value of the stress energy tensor in the relevant, 
dynamically changing, quantum state. In this way we can go beyond the standard analyses 
p8^pT| which treat the background as fixed and do not consider the non-linear quantum 
field dynamics. 

In all cases 1. - 3. , the quantum fields energy- momentum tensor is covariantly conserved 
both at the regularized as well as the renormalized levels || - [|14 



We mainly consider for the stage 3. new inflation scenarios where a scalar field evolves 
under the action of a typical symmetry breaking potential. The initial conditions will be 
taken so that the initial value of the order parameter is near the top of the potential (the 
disordered state) with essentially zero time derivative. What we find is that the existence 
of spinodal instabilities, i.e. the fact that eventually (in an expanding universe) all modes 
will act as if they have a negative mass squared, drives the quantum fluctuations to grow 
non-perturbatively large. We have the picture of an initial wave-function or density matrix 
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peaked near the unstable state and then spreading until it samples the stable vacua. Since 
these vacua are non-perturbatively far from the initial state (typically ~ to/vA, where m 
is the mass scale of the field and A the quartic self-coupling), the spinodal instabilities will 
persist until the quantum fluctuations, as encoded in the equal time two-point function 
($(x, t) 2 ), grow to 0(m 2 /X). 

This growth eventually shuts off the inflationary behavior of the scale factor as well as the 
growth of the quantum fluctuations (this last also happens in Minkowski spacetime |8|||). 

The scenario envisaged here is that of a quenched or super-cooled phase transition where 
the order parameter is zero or very small. Therefore one is led to ask: 

a) What is rolling down?. 

b) Since the quantum fluctuations are non-perturbatively large ( ~ 1/A), will not they 
modify drastically the FRW dynamics?. 

c) How can one extract (small?) metric perturbations from non-perturbatively large field 
fluctuations? 

We address the questions a)-c) as well as other issues in sec. IX. 

We choose such type of new inflationary scenario because the issue of large quantum fluc- 
tuations is particularly dramatic there. However, our methods do apply to any inflationary 
scenario as chaotic, extended and hybrid inflation. 



II. NON-EQUILIBRIUM QUANTUM FIELD THEORY, SEMICLASSICAL 

GRAVITY AND INFLATION 

We present here the framework of the non-equilibrium closed time path formalism. For 
a more complete discussion, the reader is referred to ||- ||15|| . 

The time evolution of a system is determined in the Schrodinger picture by the functional 
Liouville equation 

i^ = [H(t),P(t)i (2-1) 

where p is the density matrix and we allow for an explicitly time dependent Hamiltonian as is 
necessary to treat quantum fields in a time dependent background. Formally, the solutions to 
this equation for the time evolving density matrix are given by the time evolution operator, 
U(t,t), in the form 

p(t) = U(t,t )p(t )U- 1 (t 1 t ). (2.2) 

The quantity p(to) determines the initial condition for the evolution. We choose this initial 
condition to describe a state of local equilibrium in conformal time, which is also identified 
with the conformal adiabatic vacuum for short wavelengths. 



Given the evolution of the density matrix (|2.2|) , ensemble averages of operators are given 



by the expression (again in the Schrodinger picture) 



/nm , _ Tr[U(t ,t)OU{t,t')U(t\to)p{to)] 

{U[t)) ~ Trp(t ) ' [Z - 6) 
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where we have inserted the identity, U(t,t')U(t ,t) with t' an arbitrary time which will be 
taken to infinity. The state is first evolved forward from the initial time to to t when the 
operator is inserted. We then evolve this state forward to time t and back again to the 



initial time [B,ll 



We shall study the inflationary dynamics in a spatially flat Friedmann-Robertson- Walker 
background with scale factor a(t) and line element: 

ds 2 = dt 2 - a 2 (t) dx 2 . (2.4) 

Our Lagrangian density has the form 



(2.5) 



Our approach can be generalized to open as well as closed cosmologies. 

Our program incorporates the non-equilibrium behavior of the quantum fields involved 
in inflation into a framework where the geometry (gravity) is dynamical and is treated self 
consistently. We do this via the use of semiclassical gravity |23j where we say that the 
metric is classical and determined through the Einstein equations using the expectation 
value of the stress energy tensor (T^ u ). Such expectation value is taken in the dynamically 
determined state described by the density matrix p(t). This dynamical problem can be 
described schematically as follows: 

1. The dynamics of the scale factor a(t) is driven by the semiclassical Einstein equations 
\ + A ^ + (higher curvature) = - {T^) R . (2.6) 

OTTisrR OTlLrR 



Here Gr, Ar are the renormalized values of Newton's constant and the cosmological 
constant, respectively and G^ v is the Einstein tensor. The higher curvature terms 
must be included to absorb ultraviolet divergences. 

2. On the other hand, the density matrix pit) of the matter (that determines (T^r) 
obeys the Liouville equation ( |2.1| ) with H is the evolution Hamiltonian, which is de- 
pendent on the scale factor, a(t). 

It is this set of equations we must try to solve specifying the appropriate initial conditions. 



A. On the initial state: dynamics of phase transitions 



The situations we consider are 

1. the theory admits a symmetry breaking potential and in which the field expectation 
value starts its evolution near the unstable point. 

2. The symmetry is not broken and the field expectation value starts its evolution at a 
finite distance from the absolute minimum. 
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There is an issue as to how the field got to have an expectation value near the unstable 
point (typically at $ = 0) as well as an issue concerning the initial state of the non-zero 
momentum modes. 

Since our background is an FRW spacetime, it is spatially homogeneous and we can 
choose our state p(t) to respect this symmetry. Starting from the full quantum field t) 
we can extract its expectation value <p(t) by writing: 

$(£,*) =0(t) + *(f,t) , 0(t)=Tr[p(t)$(f,t)]^($(f,t))). (2.7) 

The quantity *$?(x, t) represents the quantum fluctuations about the zero mode (f>(t) and 
clearly satisfies (ty(x,t)) = 0. 

We need to choose a basis to represent the density matrix. A natural choice consistent 
with the translational invariance of our quantum state is that given by the Fourier modes, 
in comoving momentum space, of the quantum fluctuations ^f?(x, t): 

= J exp(-z k-x) Mt)- (2.8) 

In this language we can state our ansatz for the initial condition of the quantum state 
as follows. We take the zero mode 4>{t = 0) = O > <P(t = 0) = 0, where O w iU typically be 
very near the origin for broken symmetry and at a finite distance from it in the unbroken 
symmetry case. The initial conditions on the the nonzero modes ipk{t) will be chosen such 
that the initial density matrix p(t = 0) describes a vacuum state (i.e. an initial state in local 
thermal equilibrium at a temperature Tj = 0). There are some subtleties involved in this 



choice. First, as explained in | 12| , in order for the density matrix to commute with the initial 
Hamiltonian, we must choose the modes to be initially in the conformal adiabatic vacuum 
(these statements will be made more precise below). This choice has the added benefit of 
allowing for time independent renormalization counterterms to be used in renormalizing the 
theory. 

We are making the assumption of an initial vacuum state in order to be able to proceed 
with the calculation. It would be interesting to understand what forms of the density matrix 



can be used for other more general initial conditions [45 



The assumptions of an initial equilibrium vacuum state are essentially the same used in 



refs. ||28|| , p9[ and [|3TJ in the analysis of the quantum mechanics of inflation in a fixed de 
Sitter background. 

As discussed in the introduction, if we start from such an initial state, spinodal or para- 
metric instabilities will drive the growth of non-perturbatively large quantum fluctuations. 
In order to deal with these, we need to be able to perform calculations that take these large 
fluctuations into account. Although the quantitative features of the dynamics will depend on 
the initial state, the qualitative features associated with spinodal or parametric unstabilities 
are fairly robust for a wide choice of initial states that describe a phase transition. 



III. THE INFLATON MODEL AND THE EQUATIONS OF MOTION 

Having recognized the appearance of large quantum fluctuations driven by parametric or 
spinodal unstabilities, we need to study the dynamics within a non-perturbative framework. 
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That is, a framework allowing calculations for non-perturbatively large energy densities. 
We require that such a framework be: i) renormalizable, ii) covariant energy conserving, 
hi) numerically implementable. There are very few schemes that fulfill all of these criteria: 
the large N and the Hartree approximation ||- [14]]. Whereas the Hartree approximation is 



basically a Gaussian variational approximation |22[ that in general cannot be consistently 
improved upon, the large N approximation can be consistently implemented beyond leading 
order [^,^]. In addition, the presence of a large number of fields in most of the GUT's 
models suggest that the large N limit will be actually a realistic one. Moreover, for the case 
of broken symmetry it has the added bonus of providing many light fields (associated with 
Goldstone modes) that will permit the study of the effects of other fields which are lighter 
than the inflaton on the dynamics. Thus we will study the inflationary dynamics within 
the framework of the large N limit of a scalar theory in the vector representation of O(N) 
both for unbroken and broken symmetry. In the second case we will have a quenched phase 
transition. 



We assume that the universe is spatially flat with a metric given by eq. (|2.4|) . The matter 
action and Lagrangian density are given by eq.( |2.5|) , 



S,, 



d 4 x a 3 (t) 



■2 



-$ (x) 



1 {V${x)f 

2 a 2 (t) 



- V($(x)) 



(3.1) 




K($> T |2 + l( |2 + T-| +^^ 2 , (3-2) 

where m 2 > for unbroken symmetry and m 2 < for broken symmetry. Here lZ(t) stands 
for the scalar curvature 

m=6 (w + m, (3.3, 



a(t) a 2 (t) / 



2 



The ^-coupling of $(x) to the scalar curvature TZ(t) has been included since arises anyhow 
as a consequence of renormalization [11]. 

The gravitational sector includes the usual Einstein term in addition to a higher order 
curvature term and a cosmological constant term which are necessary to renormalize the 
theory. The action for the gravitational sector is therefore: 



S g = J d 4 x C g = J d 4 xa 3 {t) 



K(t) , « ^2 



16nG 2 



+ 77 K\t)-K 



(3.4) 



with K being the cosmological constant. In principle, we also need to include the terms 
R^Rpv and R a ^^ h 'R^^ as they are also terms of fourth order in derivatives of the metric 
(fourth adiabatic order), but the variations resulting from these terms turn out not to be 
independent of that of 1Z 2 in the flat FRW cosmology we are considering. 

The variation of the action S = S g + S m with respect to the metric gives us Einstein's 
equation 

P^ + aH^ + K g^ = -<T^> , (3.5) 

07TU 
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where is the Einstein tensor given by the variation of y/—g7Z, H^ u is the higher order 
curvature term given by the variation of yf—glZ 2 , and T^ v is the contribution from the 
matter Lagrangian. With the metric ( |2.4| ), the various components of the curvature tensors 
in terms of the scale factor are: 




-6 (ft + 3 -ft 

a 



Eventually, when we have fully renormalized the theory, we will set olr = and keep as our 
only contribution to Kr a piece related to the matter fields which we shall incorporate into 



T 



The use of semiclassical gravity (i. e. neglecting graviton loops) is justified since gravi- 
ton loops are suppressed by factors (m/Mpi) 2 where m is at the GUT scale and hence 
(m/M Pl ) 2 ~ 10" 8 . 



IV. THE LARGE N LIMIT FOR A SCALAR FIELD WITH AN ARBITRARY 

INVARIANT SELF-INTERACTION 



We present here the systematic derivation of the 1/N expansion for the scalar model 
with an arbitrary 0(iV)-invariant self-interaction and the scalar field in the vector repre- 
sentation of O(N) ||43|| - ||44|| . We consider Minkowski space-time. The generalization to the 
cosmological space-time ( |2T4] ) is discussed at the end of the section. 

The action and Lagrangian density are given by, 



S 



d x £ 



£ =2 



d$(x) ' -NV(jj $\x) 



(4.1) 



This choice of the N dependence in the interaction ensures that the large N limit exists. In 
particular, for a quartic potential we choose according to eq.(p.2|): 



A 



IN 2 



$2 



2Nm 2 
A 



The functional integral for the model takes the form 



z[J{.)\ 



expi / d 4 x \C + <S>{x).J(x 



(4.2) 



(4.3) 



where J(x) is an external source introduced to generate the correlation functions. For 
example, the connected two point correlation function is expressed as 



<T^(x) <^ k (y) >-- 



5 2 



SJi(x) 8J k (y) 



log Z[J(.)] 



(4.4) 



M-)=o 



In order to compute the large N limit it is convenient to replace the interaction term by 



the following functional integral representation [44 



S 



exp-tN J d A x V(^ 2 (x)) = J J Dz U x 5 (z(x) - ^ 2 (x)^j e~ iN I '**v[z{x)) 

= J J Dz Da e iNfd 4 x[ia(x)(z{x)~^(x))-V(z(x))} ^ 



Inserting eq.( [4.5| ) into eq.( |4.3| ), the integration over $(a;) becomes gaussian and therefore 
can be exactly computed with the result: 

)] - J J d§ Dz Da e ^ * x [?\. d$(x) Y + f <r(*)<z)-h°W v(z(x))+$(x)J(x) 
= J J Dz Da {det \d 2 + cr(.)] e iN f ^Ma^*) *0*)-V(*(*))]+5 / d 4 x d 4 y j r {x) G{x, y ,a(.)) My) ^ ^ 

Here, G(x,y,a(.)) is the inverse operator of d 2 + a(.). That is, 

~d 2 + a(x)]G(x,y,a(.)) = 5(x-y) (4.7) 

It is useful to introduce a source N K(x) for the field z(x). This will permit to generate 
the correlation functions of the composite field $ 2 (x). The generating functional now takes 
the form 

Z[J(.),K(.)]= f J Dz Da {det [d 2 + a(.)] }~ f e «r/A[^(*)^)-^W)] 

e | / rf** d 4 y Jtfa,) G?(*,»,ff(.)) J»(y)+iJV / d 4 x K(x) z(x) = f J £) z D a e iNS[z(.),a(.)] _ ( 4 _g) 

Here, 

S[z(.),ff(.)] = Jlogdet 



5 2 + a(.: 



+ / cf z 



-cr 



z(a;) — V(2(a;)) + if (a?) 



+ / rf4x ^ G (2,2/,o-(.)) ^0*0 Jifa) 



(4.9) 



Notice that the source terms in 5 are both of order one since K{x) and Ji(x) are assumed 
to be of order one. 

The functional derivatives with respect to the source K(x) at K(x) = produce the 
insertions 



1 5 



i 5K(x) 



Nz(x) = § 2 (x) 



(4.10) 



as follows from eqs. (|4.5|) and ( |4.8|) That is, K(x) is the source of the composite field $ 2 (x). 
The correlations of $ 2 (x) follow as functional derivatives of log i?[J(.), K(.)] with respect to 
K(x). For example, we have for the two points function 



<T$ 2 (x) $%) >=f 



6 2 



5K(x) 5K(y) 



logZ[J(.),K(.)) 



J(.)=,K(,)=0 



Since in eq.( |4.8| ) the iV-dependence is explicit in the exponent, we can take the large N 
limit by looking for the stationary points of the action there. 
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Extremizing S[z(.),a(.)] with respect to the field a(x) yields: 



i G(x, x, cr(.)) + z(x) 
Where we used that, 

5G(y,z,a(.)) 
Sa(x) 



1 

N 



d 4 y G(x,y,a(.)) J^y) 



1 2 



-G(y,x,a(.)) G(x,z,a(.)) . 



(4.11) 



(4.12) 



Eq. fl4.12|) can be derived as follows. Taking the functional derivative of eq. ( |4.7| ) with respect 
to <t(.) yields, 



d 2 + a(x) 



5G(x,y, <?(.)) 
5a{z) 



-S(x - z)G(x,y,a(.)) 



This is an equation for SG ^f^'^ that can be solved using the inverse operator of [d 2 + cr(x)] 
as given by eq.( |4.7| ). This gives eq. (|4.12| ). 

Extremizing S[z(.),a(.)) with respect to the field z(x) yields: 



-a(x) - V (z(x))+K(x) = . 



(4.13) 



Eqs. ( [4. 11]) -( |4. 13| ) define the saddle point a(x),z(x) as a functional of the sources 
J(.), K(.). One has thus to make the following shift of the functional integration variables 



a(x) = cr(x) + £(x) , z(x) = z(x) + w(x) 



(4.14) 



where £(x),w(x) are the new functional integration variables. 

Now we have to insert the change Q4.14 ) into eq.( |4.9|) and expand S[z(.),a(.)] in powers 
of (,(x) and w(x). The zeroth order, that is S[z(.), a(.)] yields the N = oo limit of the model. 
The quadratic part in £(x) and w(x) provides the propagators and the higher orders provide 
the vertices of the 1/N perturbation theory. 



A. The Large N Limit for the [& 2 ] 2 theory 



For simplicity, we will restrict ourselves to the $ 4 theory with potential ( |4.2j ). That is 



V(z) 



A 

8 



z + 



2m 2 V 



A 



m 



z H z + 

2 2A 



(4.15) 



In this important case we have V'(z) 
point equation ( |4.13| ), 



A z i m 



+ \ and we can eliminate z{x) using the saddle 



2 r 



z ( x ) = T a ( x ) + — m2 

A - 



(4.16) 



Equivalently, in this case we can integrate exactly over z(x) in eqs.(|4.5D, ( [4.6| ) and ( [4.8| ) since 
these functional integrals become gaussian in the field z when V(z) is given by eq. ( |4.15| ). 
The saddle point equation (|4.11| ) for u{x) becomes 
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2 r 

iG(x, x, a(.)) + — a(x) + 2K(x) - m 2 

A - 



N 



d 4 y G(x,y,a(.)) Ji(y) 



(4.17) 



where G(x, y, cr(.)) is defined by eg. (|4.71) . 

The action S at the saddle point takes then the form 



£[<)] = - log det [d 2 + a{. 



a(x) + 2K(x) - m 2 
+ 2~N J d * X d * V G(kX ' V ' a ^ ' 



(4.18) 



where we used eqs.fl4H|), (pETTD , (|4TT5|) and (FTT^) . 

Hence, from eqs.( [4.8|) and ( [4.17|) we obtain that in the N = oo limit 



lim -log^[J(.),iif(.)]=<5[a(.)]. 

N^oo iv 



(4.19) 



Furthermore, we can easily compute the two points function of the scalar field $(x) using 
eqs.([Op, QQ8D and ( CT ) with the result 



< T $ fc (j/) >= -i <5 ife G{x, y, a{.)) + O 



N 



(4.20) 



and G(x, cr(.)) was defined through eq. (f4.7[). The $(x) propagator turns to be of order 
one in the N = oo limit. 

The O(N) invariance is here explicit. We find for the expectation value of $(x) (one-point 
function) in the infinite N limit, 



6 



\ogZ[J(.),K(.)) 



5Ji(x) 

d 4 y G(x,y,a(.)) J^x) + O 



N 



6 



Ji(.)=K(.)=0 

1 

N 



SJi(x) 



S[a(.) 



Ji(.)=K(.)=0 



(4.21) 



and for the composite field $ 2 (x) 



N <*W> - — 



logZ[J(.),K(.)} 



SK(x 



S\a(.) 



.h(.)=K(.)=0 

2 



Ji(.)=K(.)=0 



A 



a(x) -m 1 +0 



(4.22) 



where we used eqs.(£OJ), ( ^T0|) and (|4l9|) . 

Once we set the external sources J;(x) and _fT(a;) equal to zero, the expectation value of 
$(a;) vanishes as it must be due to the O(N) invariance. On the contrary, N^ 1 < $ 2 (x) > 
has in general a non-zero value at zero external sources. 

It must be stressed that the present derivation applies to an arbitrary quantum state 
of the theory. Indeed, the equations simplify for the ground state due to translational 
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invariance. In such case a(x) must be a constant do and the propagator G(x,y, ao) takes 
the form, 

r d 4 k e ik -( x ~y) 

The equal-points propagator G(x, x, a ) needs an UV regulator. We obtain using a momen- 
tum cutoff A and Wick rotating k — > ik , 



G A (x,x,cr ) 



2ttH /-a k 3 dk 



(2vr) 4 7 a + k 2 (4tt) 



-A 2 + a lof 



'A 2 



1 



The A dependence can be absorbed into standard mass and coupling constant renormaliza- 
tion. Notice that p 2 = o~o is the physical (renormalized) mass squared of the fundamental 
boson in the N = oo limit as we see from eq.( f4.20|) . m 2 = m 2 B is the bare boson mass and 
A = \b the bare coupling constant. 

We find from eq.( [4.17| ) the relationship between bare and renormalized parameters, 



A 



A , 



B 



16 



^log^ 

TT Z ° 11 



B 



Ag A 2 
32vr 2 



(4.23) 



16tt 2 /i 



where A r stands for the renormalized coupling constant and we dropped all contributions of 
order A -2 and higher. 

We can now choose for simplicity Ji(x) = J(x) for % = 1, . . ., N. We proceed now to 
make a Legendre transformation such that the field expectation values <j>(x) and £(x) become 
the independent functional variables. We define following eqs.( |4.2l| ) and (|4.22j) , 



5 



p{x) 



log Z[J{.),K{.)\ = - a{x) - m 2 + 2K(x) +0[ — 



iN5K(x)'~° l " v "" WJ A 
In order to compute these functional derivatives, we used the chain rule, 



N 



(4.24) 



5 



SJ(x) 



K 



oJ(x) 



8 



K 



5cr(y) 



S 



J,K 



SJ(x) 



a,K 



The effective action functional is then given by, 

r [<K-),P(-)] = T^ \ogZ[J(.),K(.)]- Jd 4 x [J{x)4>{x)+K{x)p{x)\ 
We get for this effective action in the infinite N limit 

d 4 x 



(4.25) 



!>(.), p(.)] = - log det [d 2 + m 2 + a.) 



+ 



2A 



d 4 x d A y G(x, y, m 2 + f (.)) J(x) J(y) + O 



^x) 2 - 4^(x) 2 
1 



N 



(4.26) 



where we used eqs. (|4.18| ),( |4.19|) and (|4.24j ) and we introduced the new field variable: 
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= a(x) — m 2 = — A p(x) — 2K(x) . 



(4.27) 



Using now eqs.( |4.7|) and (|4.24 ) we can easily express J(x) in terms of <f>(x) and C,(x) as 
follows, 

1 



J{x) = d 2 + m 2 + £(x) 



(x) + 



N 



The effective action can then be written in terms of <f)(x) and £(x) with the result, 



+ 



2A 



d x 



<fx <j>{x) d 2 + m 2 + £(x) <j>(x) + O 



i{xf-AK{xf 
1 



N 



(4.28) 



The equations of motion on the fields 4>(x) and £(x) follow extremizing the effective 
action r[0(.), £(.)]. We find from eq.( |4.28| ) always for infinite N, 



d 2 + m 2 + f (z) 



(x) = 



^G(x, x, m 2 + £(.)) - \<t>{x? + ~ £(x) = . 



Or, in a more explicit form, 



d 2 + m 2 + ^</)(x) 2 - y x, m 2 + f (.)) 



(x) = 



(4.29) 



(4.30) 



d 2 + m 2 + -</>(a;) 2 - ^ x, m 2 + £(.)) 



G(x,y,m 2 + £(.)) = 6(x-y) 



As we see, this is a non-linear and non-local set of partial differential equations. 

Recall that £(x) provides the expectation value of < $ 2 (x) > as derived in eqs. (|4.22|) 
and (|427l) , 

«*) = A<«. w > + o(l 

The field <j>(x) is the expectation value of < $(x)j > for all values of % as we see from 
eqs.( P~2~TD and ( fO^) (up to C corrections). 

The ground state correspond to £(x) = <f>(x) = 0. This is the solution o" discussed above. 
We are interested on solutions with non-zero £(x) and (f>(x) describing excited states. States 
with non-zero <f>(x) will not be invariant under O(N) transformations. 



B. Invariance under spatial translations 

Let us now consider states which are invariant under spatial translations. The fields 
£(x) and <f>(x) will thus only be functions of time. This fact considerably simplifies the 
general equations (|4.30j) . We can Fourier expand G(x, y, m 2 + £(.)) as 

G(x,y,m 2 + £(.)) = / e*<*-« g k (t,t') 
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where the times t and t' are associated with x and y respectively and gk{t,t') obeys the 
equation 



dt 2 



+ k 2 + m 2 + f (t) 



g h (t,t') = 5(t-t') . 



(4.31) 



This one-dimensional Green function can be expressed in terms of solutions of the ho- 
mogeneous equation 



dt 2 



+ k 2 + m 2 + £(t) 



fk{t) = • 



(4.32) 



One gets, 



9k(t,t') 



#(*<)#(*>) 



where t< = min(i, £'), t > = max(t,t') and f k (t) and f^(t) are independent solutions of 
eq. (|4.32|) . W k stands for the Wronskian between these solutions: 



W k = W{f<(t)J>(t)}. 

For causal bondary conditions one has f k {t) = [f k (t)}" 
Before renormalization, eqs. (|4.30|) take the form 



(4.33) 



d 2 2 , ^jl/j.\2 



„ + m 2 + -6(t) 
dt 2 2 rw 2 



where G(x, x, m 2 + £(.)) = J 



G(x,x,m 2 + £(.)) 
d 3 k \f k (t)\ 



<f>(t) = 



(2tt)3 W k 



(4.34) 



and the mode functions f k (t) obey, 



1_ + k 2 + m 2 + ^0( t )2 _ ^ G(Xj Xj m 2 + £Q) 



/*(*) = o 



(4.35) 



We obtained an infinite set of coupled ordinary differential equations on <p(t), f k (t) for 
< k < oo. Notice that eqs.( [4.34"D -( |4.35| ) are local on time. That is they involve the 
unknown functions <p(t), f k (t) always at time t. 

All physical quantities can be computed in terms of the mode functions f k {t), < k < oo 
and the order parameter 4>(t). We have described the large N approximation. Other non- 
perturbative approximations are the Hartree approximations and the self-consistent one-loop 
approximation. They were considered in ref. 0. 



In a cosmological spacetime (|2.4j) the large N evolution equations take the form 

4>(t) = o 



&% + 3H{t) Jt + mt) + m2 + + ^ 2{t)) 



d 2 nrrf \d 



k2 +mt)+m 2 + ^(t) 2 + ^ 2 (t)} 



fk(t) = 



(4.36) 
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where 



H(t) = 



ix(t) 
a(tj 



and we introduced the notation 

(< n *(t)) = -iG(x,x,m 2 + £(.)) 



d 3 k \f k (t)f 



(2vr) 3 iW k 



(4.37) 



C. The Large N Limit in Cosmological Spacetimes 

Let us connect the functional formalism in previous sections with the operatorial ap- 
proach. Let us call '1' the direction in internal space of the expectation value of $(x). We 
write, 

$(x,t) = (a(x,t),Tf(x,t)), 

with 7? an (N — l)-plet, and 

a(x,t) = VN<P(t) + X (x,t) ; (a(x, t)) = VN<j>(t) ; ( X (x,t)) = 0. (4.38) 

We can now write the operator tt(x, t) in terms of creation and annihilation operators and 
mode functions that obey the Heisenberg equations of motion 



tt(x, t) 



d 3 k 



a k f k (t) e ik - s + 4 f* k (t) e~ lk 



ik-x 



(4.39) 



where a k and ot obey canonical commutation rules. It is easy to derive the two points 



function eq.(^20|) from eqs.(|L3g), ( ^33|) and Q439|) . 

We see that since there are N — 1 'pion' fields, contributions from the field x can be 
neglected in the N — > oo limit as they are of order 1/N with respect those of 7r and <j). 

The equations of motion ( fL36|) can be written as 



where 



4>{t) + 3H(t) <j>{t) + M 2 {t) 0(t) = 



M\t) = m 2 + £ Tl{t) + ^0 2 (t) + ^(vr 2 (t)> . 



(4.40) 



(4.41) 



plays the role off time-dependent effective mass. 

An important point to note in the large N equations of motion is that the form of the 
equation for the zero mode ( 4.40|) is the same as for the k = mode function (|4.40 ). It is this 
property that allows solutions of these equations in a symmetry broken scenario to satisfy 
Goldstone's theorem in out of equilibrium situations both in Minkowski and cosmological 
spacetimes 0,00,111 . 
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In this leading order in 1/N the theory becomes Gaussian, but with the self-consistency 
condition (|4.37|) . 

The initial conditions on the modes /&(£) must now be determined. At this stage it 
proves illuminating to pass to conformal time variables in terms of the conformally rescaled 
fields 



sec 



12 ] and section V for a discussion) in which the mode functions obey an equation 
which is very similar to that of harmonic oscillators with time dependent frequencies in 
Minkowski space-time. It has been realized that different initial conditions on the mode 
functions lead to different renormalization counterterms |12[]; in particular imposing initial 
conditions in comoving time leads to counterterms that depend on these initial conditions. 
Thus we chose to impose initial conditions in conformal time in terms of the conformally 
rescaled mode functions leading to the following choice: 



with initial conditions in comoving time, 

1 



fk(k 



with 



fk(to) 



d(t c 



k 2 + M 2 (t )- 



a(to) 
K(t Q ) 



fk(tc 



6 



(4.42) 



(4.43) 



We thus find from eq. (|4.33|) that the Wronskian takes the value Wk = —2i and the quantum 
fluctuations of the inflaton ( f4.37| ) take the form 



d 3 k 

2(2vr) ; 



\fk(t)t 



(4.44) 



For convenience, we have set a(io) = 1 in eq. ([4.43 ). At this point we recognize that when 
■M 2 (t ) — lZ(t )/6 < the above initial condition must be modified to avoid imaginary 
frequencies, which are the signal of instabilities for long wavelength modes in the broken 
symmetry case. Thus we define the initial frequencies that determine the initial conditions 
( TO as 



M 2 (t c 



k 2 + 

n 2 k = k 2 + M 2 (t c 



K(t ) 



6 

K(t ) 
6 



for k 2 < 



M 2 (t ) - 



K(t ) 



for k 2 > 



M 2 (t ) 



K{t c 



6 



(4.45) 
(4.46) 



In the unbroken symmetry case (m > ) we use eq.( |4.46| ) for all k 



As an alternative we have also used initial conditions which smoothly interpolate from 
positive frequencies for the unstable modes to the adiabatic vacuum initial conditions defined 

While the alternative choices of initial conditions 
(a few percent in quantities which 



by ( [4.42|) -( ^43|) for the high k modes, 
result in small quantitative differences in the results 
depend strongly on these low-k modes), all of the qualitative features we will examine are 
independent of this choice. 
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In the large N limit we find the energy density and pressure density to be given by [p~3| . p~2 



e — 3p 

N 
1 



, , 2m 4 
+ 2m z 4> z + - <f + — - 
2 A 



(7T 2 ) 



+ ^((Vvr) 2 ) + 2m 2 (vr 2 ) - £ ^ (tt 2 ) + ^ [20 2 (vr 2 ) + (vr 2 ) 2 ] + 6 £ (<tt tt> + (tV 2 ) + 3 £ (tttt)) 



where (7r 2 ) is given by equation ( [4.44 ) and we have defined the following integrals: 



((Vtt) 2 ) 



d 3 A; 
2(2^ 



k 2 \h{t)\ 2 , (tt 2 ) 



2(2tt) 



|A(*)I S 



(4.47) 



The composite operators (ir-h) and (tttt) are symmetrized by removing a normal ordering 
constant, (tttt) may be rewritten using the equation of motion ( J4.4U ) : 



(rft) = - 3 * ( tttt) - - M 2 (t) (tt 2 ) 



(4.46 



It is straightforward to show that the bare energy is covariantly conserved by using the 
equations of motion for the zero mode and the mode functions. 



V. RENORMALIZATION, CONFORMAL TIME AND INITIAL CONDITIONS 



Renormalization is a very subtle but important issue in gravitational backgrounds 
The fluctuation contribution (tc 2 (x, t)), the energy, and the pressure all need to be renormal- 
ized. The renormalization aspects in curved space times have been discussed at length in the 
literature |23| and have been extended to the large N self-consistent approximations for the 
non-equilibrium backreaction problem in |K^ . |l^ . |T2| . |33|j3^] . More recently, a consistent and 
covariant regularization scheme that can be implemented numerically has been proposed 



24 



The ultraviolet divergences can be seen in the present framework in the fc-integrals over 



the modes in eqs.( [4.44"D and ( [4.47] ). To analyze these divergences it is convenient to change 



t fit' 

T = ?o+ I -77^ I T(t = to) = T , (5.1) 



variables to conformal time T defined as 

ft dt' 
ho a(t' 

The metric becomes then 

ds 2 = C 2 (T) (dT 2 - dx 2 ) , (5.2) 

where C(T) = a(t(T)) ; C(T ) = 1 stands for the scale factor in conformal time. 

The issue of renormalization and initial conditions is best analyzed in conformal time 
which is a natural framework for adiabatic renormalization and regularization. 
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Under a conformal rescaling of the field 



$(f, t) 



C(T) 



(5.3) 



the action for a scalar field (with the obvious generalization to iV components) becomes, 
after an integration by parts and dropping a surface term 



S 



d 3 x dT 



\{xf-\iyxf-v{x) 



with 



V(x) = C\T) V 



X 



C{T), 



C\T) g x 2 



(5.4) 



(5.5) 



where TZ = 6 C"(T)/C 3 (T) is the Ricci scalar, and primes stand for derivatives with respect 
to conformal time T. 

The conformal time Hamiltonian operator, which is the generator of translations in T, 
is given by 



Hr = /^{^ + i(V X ) 2 + V(x) 



(5.6) 



with Il x being the canonical momentum conjugate to x-, n x = x' ■ Separating the zero mode 
of the field x 



X{x, T) = XoCO + x{x, T), 



(5.7) 



and in the large approximation we find that the Hamiltonian becomes linear plus quadratic 
in the fluctuations, and similar to a Minkowski space-time Hamiltonian with a T-dependent 
effective mass term given by 



B(T) = C 2 (T) 



™ 2 + tt-\)K + \xl(T) + \(x 2 ) 



(5.8) 



Notice that this B(t) naturally appears in the WKB expansion of the mode functions both 
in cosmic and conformal time [see eq.( 5.24 )]. 



We can now follow the steps and use the results of reference [T1J for the conformal time 
evolution of the density matrix by setting a(t) = 1 in the proper equations of that reference 
and replacing the frequencies by 

4(T) = k 2 +B(T), (5.9) 

and the expectation value in eq. (|5.8j) is obtained in this T evolved density matrix. The time 
evolution of the kernels in the density matrix (see JOJ) is determined by the mode functions 
that obey 



' d 2 
dT 2 



+ k 2 + B{T) 



F k (T) = . 



(5.10) 
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The Wronskian of these mode functions 



W(F, F* 



F' k F* k - F k F* 



(5.11) 



is a constant. It is natural to impose initial conditions such that at the initial T the density 
matrix describes a situation of local thermodynamic equilibrium and therefore commutes 
with the conformal time Hamiltonian at the initial time. This implies that the initial con- 
ditions of the mode functions F k (T) be chosen to be (see [|TT|]) 



F k (T ) 



u k (T ) 



; F k( T o) = -iJuj k (T ) . 



With such initial conditions, the Wronskian ( f5.11| ) takes the value 

W(F, F*) = -2% . 



(5.12) 



(5.13) 



These initial conditions correspond to the choice of mode functions which coincide with the 
first order adiabatic modes and those of the Bunch-Davies vacuum for large momentum 
To see this clearly, we write the solution of eq. ( |5.10|) in the form, 



D k {T) = jl Rk{r ' )dT \ 
with the function R k (T) obeying the Riccati equation 

R' k + R 2 k + k 2 + B(T) = . 
This equation posses the solution 



R k (T) = -ik + Ro, k (T) - i 



Ri,k(T) . #2,fc(2~) .R3,k(T) R^ k {T) 



k k 2 k 3 

and its complex conjugate. We find for the coefficients: 



+ 



k 4 



+ o 



(5.14) 



(5.15) 



(5.16) 



Ro,k — ; R\ tk — - B(T) ; R 2}k 



2 R ' l > k 



^3,k — p (R'2,k ~ R\,k) i -^4,fe — -p (^3,fe + 2 -Rl,fc-R2,fc 



(5.17) 



The solutions F k (T) obeying the boundary conditions (|5.12 ) are obtained as linear combi- 
nations of this WKB solution and its complex conjugate 



F k (T) 



1 



2Ju k (T ) 



[(i + 7 )L> fe (r) + (i- 7 )^(r)] , 



(5.18) 



where the coefficient 7 is obtained from the initial conditions. It is straightforward to find 
that the real and imaginary parts are given by 



lR = 1 + 0(l/k 4 ) ; 7/ = 0(l/k 3 ) . 



(5.19) 



Therefore the large-fc mode functions satisfy the adiabatic vacuum initial conditions ||23| 
This, in fact, is the rationale for the choice of the initial conditions ( 5.12 ). 
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Following the analysis presented in [|11] we find, in conformal time that 

d 3 k 



2(2n) 3 



\Fh(T)\' 



(5.20) 



The Heisenberg field operators x(af, T) and their canonical momenta U x (x,T) can be ex- 
panded as: 



x{x, T) 



U x (x,T) 



d 3 k 



x/2(2vr) 3 / 2 
d 3 k 



a k F k (T) + al k F*{T) 
~a k F^T)+al k FZ'(T) 



ik-x 



ik-x 



(5.21) 
(5.22) 



with the time independent creation and annihilation operators a k and a k obeying canonical 
commutation relations. Since the fluctuation fields in comoving and conformal time are 
related by the conformal rescaling (|5.3j ), it is straightforward to see that the mode functions 
in comoving time are related to those in conformal time simply as 



/*(*) 



F k {T) 
C(T)- 



(5.23) 



Therefore the initial conditions ( |5.12j ) on the conformal time mode functions imply the initial 
conditions for the mode functions in comoving time are given by eq. ( 14.421) . 

For renormalization purposes we need the \axge-k behavior of \ fk{t)\ 2 , \fk{t)\ 2 , which are 
determined by the large-fc behavior of the conformal time mode functions and its derivative. 
These are given by appropriately adapting the Minkowski formulas |§ 



\FJT) 



|f*COl a 



k 2 k 2 



fc 5 



k + 



B{T) 1 
2k 8 k 3 



3B 2 (T) + B"(T) 
B 2 (T) + B"(T) 



0, F 



(5.24) 



We note that the large k behavior of the mode functions to the order needed to renormalize 
the quadratic and logarithmic divergences is insensitive to the initial conditions. This is not 
the case when the initial conditions are imposed as described in [II,I3|. Thus the merit in 



considering the initial conditions in conformal time [12 



There is an important physical consequence of the choice ( 5.12j ) of initial conditions, 
which is revealed by analyzing the evolution of the density matrix. 

In the large iV or Hartree (also to one-loop) approximation, the density matrix is Gaus- 
sian, and defined by a normalization factor, a complex covariance that determines the diag- 
onal matrix elements and a real covariance that determines the mixing in the Schrodinger 
representation as discussed in reference 0] (and references therein). 

In conformal time quantization and in the Schrodinger representation in which the field 
X is diagonal the conformal time evolution of the density matrix is via the conformal time 
Hamiltonian (|5.6| ). The evolution equations for the covariances is obtained from those given 
in reference QTTJ by setting a(t) = 1 and using the frequencies u) k (T) = k 2 + A4 2 (T). In 
particular, by setting the covariance of the diagonal elements (given by equation (2.20) in 



|TT| ; see also equation (2.44) of jlTI), 
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A k {T) 



. F' k *(T) 
'WO 



(5.25) 



we find that with the initial conditions ( |5.12| ), the conformal time density matrix is that 
of local equilibrium at % in the sense that it commutes with the conformal time Hamilto- 
nian. However, it is straightforward to see, that the comoving time density matrix does not 
commute with the comoving time Hamiltonian at the initial time to- 

An important corollary of this analysis and comparison with other initial conditions 
used in comoving time is that assuming initial conditions of local equilibrium in comoving 
time leads to divergences that depend on the initial condition as discussed at length in 
This dependence of the renormalization counterterms on the initial condition was 
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also realized in ref. |2T| within the context of the CTP formulation. Imposing the initial 
conditions corresponding to local thermal equilibrium in conformal time, we see that: i) 
the renormalization counterterms do not depend on the initial conditions and ii) the mode 
functions are identified with those corresponding to the adiabatic vacuum for large momenta. 
This is why we prefer the initial conditions ( |5.12j ). 

For our main analysis we choose this initial temperature to be zero so that the resulting 
density matrix describes a pure state, which for the large momentum modes coincides with 
the conformal adiabatic vacuum. Such zero temperature choice seems appropriate after the 
exponential inflation of the universe. 

Particle Number: 

We write the Fourier components of the field \ and its canonical momentum Il x given 
by (^ -(£23) as: 



V2 



a k F k (T)+al k F*(T) 



U Xik (T) = -j=[a k F^T) + a[ 



(5.26) 
(5.27) 



These (conformal time) Heisenberg operators can be written equivalently in terms of the T 
dependent creation and annihilation operators 



n*,*(T) = -i 



~a k (T) e - l ^ ro)T + 4(T) e iuJk{To)r 
~a k (T) e - iuJk{To)T - 4(T) e luJk{To)T 



(5.28) 
(5.29) 



The operators a k {T) ;a k (T) are related by a Bogoliubov transformation. The number of 
particles referred to the initial Fock vacuum of the modes F k , is given by 



N k (T) = (al(T)~a k (T)) = - 



F k (%) 



+ 



F' k {T) 



F k (%) 



(5.30) 



or alternatively, in terms of the comoving mode functions f k (t) = F k (T)/C(T) we find 



N k (t) 



a 2 {t) 



f k (t ) 



+ 



v k (h 



f k (t) + H(t)f k (t) 
fk(to) 



1 

2 



(5.31) 
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Using the large fc-expansion of the conformal mode functions given by eqs. ( |5.24j) we find 
the large-/c behavior of the particle number to be Nk fc =°° 0{l/k 4 ), and the total number of 
particles (with reference to the initial state at T ) is therefore finite. 
Renormalization. 

We make our subtractions using an ultraviolet cutoff, Aa(t), constant in physical co- 
ordinates. This guarantees that the counterterms will be time independent [|TJ|]. The 
renormalization then proceeds much in the same manner as in reference ||11|| ; the quadratic 
divergences renormalize the mass and the logarithmic terms renormalize the quartic coupling 
and the coupling to the Ricci scalar. In addition, there is a quartic divergence which renor- 
malizes the cosmological constant while the leading renormalizations of Newton's constant 
and the higher order curvature coupling are quadratic and logarithmic respectively. The 
renormalization conditions on the mass, coupling to the Ricci scalar and coupling constant 
are obtained from the requirement that the frequencies that appear in the mode equations 
are finite Ol, i.e: 



A 



A 



A 



A 



ml + £ B K(t) + ^\t) + '-f(n 2 (t)) B = m 2 R + £ R K(t) + f<P 2 (t) + 'f(n\t)) R , (5.32) 

while the renormalizations of Newton's constant, the higher order curvature coupling, and 
the cosmological constant are given by the condition of finiteness of the semi-classical 
Einstein-Friedmann equation: 



r " + a B H° + K B g° + (T °) B = + a R H° + K R g° + (T°) R . 



8ttG b u " u n u/ 8hG r 

Finally, we arrive at the following set of renormalizations [ |T4"| , |T2" 
1 



1 



8ttNG r 8ttNG b (4tt) 2 
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(5.33) 



Here, k is the renormalization point. As expected, the logarithmic terms are consistent with 
the renormalizations found using dimensional regularization [24"|,54"||. Again, we set a R = 
and choose the renormalized cosmological constant such that the vacuum energy is zero in 
the true vacuum. We emphasize that while the regulator we have chosen does not respect 
the covariance of the theory, the renormalized energy momentum tensor defined in this way 
nevertheless retains the property of covariant conservation in the limit when the cutoff is 
taken to infinity. 

The logarithmic subtractions can be neglected because of the coupling A < 1CT 12 . Using 
the Planck scale as the cutoff and the inflaton mass m R as a renormalization point, these 
terms are of order A ln[M p i/m R ] < 10~ 10 , for Mpi > m > 10 9 GeV . An equivalent statement 
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is that for these values of the coupling and inflaton masses, the Landau pole is well beyond the 
physical cutoff M p i. Our relative error in the numerical analysis is of order 10" 8 , therefore our 
numerical study is insensitive to the logarithmic corrections. Though these corrections are 
fundamentally important, numerically they can be neglected. Therefore, in the numerical 
computations that follow, we will neglect logarithmic renormalization and subtract only 
quartic and quadratic divergences in the energy and pressure, and quadratic divergences in 
the fluctuation contribution. 



A. Renormalized Equations of Motion for Dynamical Evolution in the Large N limit 

It is convenient to introduce the following dimensionless quantities and definitions, 

r = m R t ; h(r) = ; q= u q = ; g = — , (5.34) 

m R m R m R 8n z 
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(n 2 (t)) R ; f q (r) = y/^ f k (t) . (5.35) 



Choosing £ R = (minimal coupling) and the renormalization point k = \m R \ and setting 
a(r ) = 1, the equations of motion become for unbroken symmetry: 
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We find for broken symmetry, 



r)(r) = , 

f q (r) = , 



(5.37) 
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(5.39) 



both for unbroken and broken symmetry. 

The initial conditions for t](t) will be specified later. An important point to notice is 
that the equation of motion for the q = mode coincides with that of the zero mode ( p.37| ). 
Furthermore, for r/(r — > oo) 7^ 0, a stationary (equilibrium) solution of the eq.( |5.37| ) is 
obtained for broken symmetry when the sum rule PJ5| jr3]JT2] | 



1 +r/ 2 (oo) + #£(00) = 



(5.40) 



is fulfilled. 

This sum rule is nothing but a proof that Goldstone's theorem holds here out of thermal 
equilibrium. In addition, it is a result of the fact that the large N approximation satisfies 
the Ward identities associated with the O(N) symmetry, since the term —1 + rf{r) + gE(r) 
is seen to be the effective mass of the modes transverse to the symmetry breaking direction, 
i.e. the Goldstone modes in the broken symmetry phase. 

The renormalized dimensionless evolution equations in the Hartree approximation are 
very similar to eqs. (|5.36|) - (|5.38|) . They can be obtained just dividing by three the r/ 2 (r) term 
in the zero mode equation ||. 

In terms of the zero mode 7/(t) and the quantum mode function given by eq.( |5.37| ) we find 
that the Friedmann equation for the dynamics of the scale factor in dimensionless variables 
is given by 



h 2 {r)=Ahl e R (r) 
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and the renormalized energy and pressure are given by: 
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where the subtractions <S^ and are given by 
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The renormalized energy and pressure are covariantly conserved: 
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e B (r) + 3/i(r)(p + e)fl(r) = 0. (5.43) 

From the evolution of the mode functions that determine the quantum fluctuations, we can 
study the growth of correlated domains with the equal time correlation function, 

S(x,t) = (vr(f,t)vr(0,t)) = J e** |/ fc (t)| 2 , (5.44) 

which can be written in terms of the power spectrum of quantum fluctuations, |/ 9 (t)| 2 . It 
is convenient to define the dimensionless correlation function, 

S(\x\,t) 1 r°° 

s{p, T ) = — ^— = ir^T / 1 d( l si %p] \ fqi. r )\ 2 ! P = m n\x\ ■ ( 5 - 45 ) 

m R ^TT z p Jo 

We now have all the ingredients to study the particular cases of interest. 



VI. SCALAR FIELD DYNAMICS IN A FIXED FRW BACKGROUND 



We consider in this section the evolution of scalar fields in radiation or matter dominated 
FRW cosmologies |12| . The case for de Sitter expansion will be discussed in sec. IX |13| . 

We write the scale factor as a(r) = (r/ro) n with n = 1/2 and n = 2/3 corresponding to 
radiation and matter dominated backgrounds, respectively. The value of To determines the 
initial Hubble constant since 

a(r ) _ n_ 
a(r Q ) r ' 



h(r ) 



We now solve the system of equations ( 4.36|) in the large N limit. We begin by presenting 
an early time analysis of the slow roll scenario. We then undertake a thorough numerical 
investigation of various cases of interest. For the symmetry broken case, we also provide an 
investigation of the late time behavior of the zero mode and the quantum fluctuations. We 
use the dimensionless variables (|5.34| )- ([5735| ). 

We will assume minimal coupling to the curvature, £ r = 0. In the cases of interest, 
1Z <^ /i 2 , so that finite £ r has little effect. 



A. Early Time Solutions for Slow Roll 

For early times in a slow roll scenario [m 2 = — fi 2 , r](r ) 1], we can neglect in eqs.( ^.38|) 
both the quadratic and cubic terms in r/(r) as well as the quantum fluctuations (7r 2 (r)) r 
[recall that (7r 2 (ro)) r = 0]. Thus, the differential equations for the zero mode and the mode 
functions Q4.36| ) become linear equations. In terms of the scaled variables introduced above, 
with a(r) = r n (n = 2/3 for a matter dominated cosmology while n = 1/2 for a radiation 
dominated cosmology) we have: 

3?7< 

V{r) + — j)(t)-»j(t) = 0, (6.1) 
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The solutions to the zero mode equation ( |6.1|) are 

r 1 {r) = cr- u I u {T) + dT- u K u {T) 



(6.3) 



where v = (3n — l)/2, and I v {t) and K v (t) are modified Bessel functions. The coefficients, 
c and d, are determined by the initial conditions on 77. For f](r ) = i] and ?)(t ) = 0, we 
have: 
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Taking the asymptotic forms of the modified Bessel functions, we find that for intermediate 
times f]{r) grows as 



v{t) ^ _JL r -3n/2 , 
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(6.4) 



We see that 77 (r) grows very quickly in time, and the approximations ( |6.1|) and ( |6.2|) will 
quickly break down. For the case shown in fig.l (with n = 2/3, ^(to) = 10^ 7 , and tj{to) = 0), 
we find that this approximation is valid up to r — tq ~ 10. 

The equations for the mode functions ( |6.2|) can be solved in closed form for the modes 
in the case of a radiation dominated cosmology with n = 1/2. The solutions are 
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Here, [/(•) and M(-) are confluent hypergeometric functions |25[ (in another common nota- 
tion, M(-) = iFi(-)), and the c k and (i^ are coefficients determined by the initial conditions 
( fl.42| ) on the modes. The solutions can also be written in terms of parabolic cylinder func- 
tions. 

For large r we have the asymptotic form 
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(6.6) 



Again, these expressions only apply for intermediate times before the nonlinearities have 
grown significantly. 



B. Numerical Analysis 

We now present the numerical analysis of the dynamical evolution of scalar fields in time 
dependent, matter and radiation dominated cosmo logical backgrounds. We use initial values 
of the Hubble constant such that h(ro) > 0.1. For expansion rates much less than this value 
the evolution will look similar to Minkowski space, which has been studied in great detail 
elsewhere lOf] . As will be seen, the equation of state found numerically is, in the majority 



of cases, that of cold matter. We therefore use matter dominated expansion for the evolution 
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in much of the analysis that follows. The evolution in radiation dominated universes remains 
largely unchanged, although there is greater initial growth of quantum fluctuations due to 
the scale factor growing more slowly in time. Using the large N approximation to study 
theories with continuous and discrete symmetries respectively, we treat three important 
cases. They are 1) m 2 < 0, r/(r ) < 1; 2) m 2 < 0, r/(r ) > 1; 3) m 2 > 0, 77(70) > I. 

In presenting the figures, we have shifted the origin of time such that r — > r' = r — tq. 
This places the initial time, To, at the origin. In these shifted coordinates, the scale factor 
is given by 

a{t) = 

V r 

where, once again, n = 2/3 and n = 1/2 in matter and radiation dominated backgrounds 
respectively, and the value of r is determined by the Hubble constant at the initial time: 

h(r ) = -. 

Case 1: m 2 < 0, tj(tq) <C 1. This is the case of an early universe phase transition in 
which there is little or no biasing in the initial configuration (by biasing we mean that the 
initial conditions break the 77 — > —77 symmetry). The transition occurs from an initial tem- 
perature above the critical temperature, T > T c , which is quenched at to to the temperature 
Tf <C T c . This change in temperature due to the rapid expansion of the universe is modeled 
here by an instantaneous change in the mass from an initial value m 2 = T 2 /T 2 — 1 to a final 
value m 2 = —1. We will use the value m 2 = 1 in what follows. This quench approximation 
is necessary since the low momentum frequencies (|4.43|) appearing in our initial conditions 
( |4.42j ) are complex for negative mass squared and small t]{to). An alternative choice is to 
use initial frequencies given by 



wjfe(r ) 



k 2 + M 2 (t ) tanh 



' fc 2 + A4 2 (r ) Nll/2 
, \M 2 (t )\ , 



These frequencies have the attractive feature that they match the conformal adiabatic fre- 
quencies given by eq.( 4.43| ) for large values of k while remaining positive for small k. We 



find that such a choice of initial conditions changes the quantitative value of the particle 
number by a few percent, but leaves the qualitative results unchanged. 

We plot the the zero mode T](t), the equal time correlator grS(r), the total number of 
produced particles gN(r) (see sec. VI for a discussion of our definition of particles), the 
number of particles gN^r) as a function of wavenumber for both intermediate and late 
times, and the ratio of the pressure and energy densities p(r)/e(r) (giving the equation of 
state). 

Figs, la-e shows these quantities in the large iV approximation for a matter dominated 
cosmology with an initial condition on the zero mode given by 77(70) = 10 -7 , 77(70) = 
and for an initial expansion rate of H(tq) = 0.1. This choice for the initial value of r] stems 
from the fact that the quantum fluctuations only have time to grow significantly for initial 
values satisfying T](tq) ^fg\ for values 77(70) ■Jg the evolution is essentially classical. 
This result is clear from the intermediate time dependence of the zero mode and the low 
momentum mode functions given by the expressions (|6.4j) and (|6.6|) respectively. 
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FIGURES 



Figure la Figure lb 




100 200 300 400 t 



FIG. 1. Symmetry broken, slow roll, large N, matter dominated evolution of (a) the zero 
mode rj(r) vs. r, (b) the quantum fluctuation operator <?£(t) vs. t, (c) the number of particles 
gN(r) vs. r, (d) the particle distribution gNk(r) vs. k at r = 149.1 (dashed line) and r = 398.2 
(solid line), and (e) the ratio of the pressure and energy density p(r)/e(r) vs. t for the values, 
rj(r ) = 10- 7 , 77(70) = 0, g = 1(T 12 , h(r ) = 0.1. 

After the initial growth of the fluctuation gS(r) (fig. lb) we see that the zero mode 
(fig. la) approaches the value given by the minimum of the tree level potential, 77 = 1, while 
gE(r) decays for late times as 

c c 
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For these late times, the Ward identity corresponding to the O(N) symmetry of the field 
theory is satisfied, enforcing the condition 



-l +v 2 (r)+g^(r) = 0. (6.7) 
Hence, the zero mode approaches the classical minimum as 

C 



7] (t) ~ 1 - 



a?{r) 



Figure lc depicts the number of particles produced. After an initial burst of particle 
production, the number of particles settles down to a relatively constant value. Notice that 
the number of particles produced is approximately of order 1/g. In fig. Id, we show the 
number of particles as a function of the wavenumber, k. For intermediate times we see the 
simple structure depicted by the dashed line in the figure, while for late times this quantity 
becomes concentrated more at low values of the momentum k. 

Finally, fig.le shows that the field begins with a de Sitter equation of state p = — e but 
evolves quickly to a state dominated by ordinary matter, with an equation of state (averaged 
over the oscillation timescale) p = 0. This last result is a bit surprising as one expects from 
the condition ( |6.7|) that the particles produced in the final state are massless Goldstone 
bosons which should have the equation of state of radiation. However, as shown in fig. Id, 
the produced particles are of low momentum, g C 1, and while the effective mass of the 
particles is zero to very high accuracy when averaged over the oscillation timescale, the 
effective mass makes small oscillations about zero so that the dispersion relation for these 
particles differs from that of radiation. In addition, since the produced particles have little 
energy, the contribution to the energy density from the zero mode, which contributes to a 
cold matter equation of state, remains significant. 

Finally, we show the special case in which there is no initial biasing in the field, i](t ) = 
0, ?)(ro) = 0, and /i(ro) = 0.1 in figs. 2a-d. The zero mode remains zero for all time, so that 
the quantity g£(r) (fig. 2a) satisfies the sum rule ( |6.7D by reaching the value one without 
decaying for late times. Notice that many more particles are produced in this case (fig 2b); 
the growth of the particle number for late times is due to the expansion of the universe. The 
particle distribution (fig. 2c) is similar to that of the slow roll case in fig.l. The equation of 
state (fig. 2d) is likewise similar. 

In each of these cases of slow roll dynamics, increasing the Hubble constant has the 
effect of slowing the growth of both rj and g^{t). The equation of state will be that of a de 
Sitter universe for a longer period before moving to a matter dominated equation of state. 
Otherwise, the dynamics is much the same as in figs. 1-3. 

Case 2: m 2 < 0, ^(to) ^> 1. We now examine the case of a chaotic inflationary scenario 
with a symmetry broken potential. In chaotic inflation, the zero mode begins with a value 
rj(r) ^> 1. During the de Sitter phase, ft > 1, and the field initially evolves classically, 
dominated by the first order derivative term appearing in the zero mode equation [see 
eq . (|4.36|) 1 . Eventually, the zero mode rolls down the potential, ending the de Sitter phase 



and beginning the FRW phase. We consider the field dynamics in the FRW universe after 
the end of inflation. We thus take the initial temperature to be zero, T = 0. 
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Figure 2a Figure 2b 
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FIG. 2. Symmetry broken, r)(j) = 0, matter dominated evolution of (a) the quantum fluctua- 
tion operator gS(r) vs. r, (b) the number of particles gN{r) vs. r, (c) the particle distribution 
gNk(r) vs. k at r = 150.1 (dashed line) and r = 397.1 (solid line), and (d) the ratio of the pressure 
and energy density p(r)/e(r) vs. r for the values ??(ro) = 0, f)(ro) = 0, g = 10 -12 , /i(to) = 0.1. 

Figure 4 shows our results for the quantities, rj{r), gX(r), gN(r), gN k {j) and p{r)/e{r) 
for the evolution in the large iV approximation within a radiation dominated gravitational 
background with h{ro) = 0.1. The initial condition on the zero mode is chosen to have the 
representative value //(to) = 4 with i)(tq) = 0. Initial values of the zero mode much smaller 
than this will not produce significant growth of quantum fluctuations; initial values larger 
than this produces qualitatively similar results, although the resulting number of particles 
will be greater and the time it takes for the zero mode to settle into its asymptotic state 
will be longer. 

We see from fig. 3a that the zero mode oscillates rapidly, while the amplitude of the 
oscillation decreases due to the expansion of the universe. 

This oscillation induces particle production through the process of parametric amplifi- 
cation (fig. 3c) and causes the fluctuation (?X(t) to grow (fig.3b). Eventually, the zero mode 
loses enough energy that it is restricted to one of the two minima of the tree level effective 
potential. The subsequent evolution closely follows that of Case 1 above with g£(r) decay- 
ing in time as l/a 2 (r) ~ 1/r with t}(t) given by the sum rule (|6.7|). The spectrum (fig. 3d) 
indicates a single unstable band of particle production dominated by the modes k — 1/2 
to about k = 3 for late times. The structure within this band becomes more complex with 
time and shifts somewhat toward lower momentum modes. 



30 




FIG. 3. Symmetry broken, chaotic, large N, radiation dominated evolution of (a) the zero 
mode ry(r) vs. r, (b) the quantum fluctuation operator gS(r) vs. r, (c) the number of particles 
gN(r) vs. r, (d) the particle distribution gNk(r) vs. k at r = 76.4 (dashed line) and r = 392.8 
(solid line), and (e) the ratio of the pressure and energy density p(r)/e(r) vs. r for the values 
77(70) = 4, 77(70) = 0, g = 10- 12 , fc(7o) = 0.1. 

Such a shift is also observed in Minkowski spacetimes |||9j . Figure 4e shows the equation 
of state which we see to be somewhere between the relations for matter and radiation 
for times out as far as t = 400, but slowly moving to a matter equation of state. Since 
matter redshifts as l/a 3 (t) while radiation redshifts as l/a 4 (rj), the equation of state should 
eventually become matter dominated. Given the equation of state indicated by fig.3e, we 
estimate that this occurs for times of order t = 10 4 . The reason the equation of state in this 
case differs from that of cold matter as was seen in figs. 1-3 is that the particle distribution 
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FIG. 4. Symmetry broken, chaotic, large N, matter dominated evolution of (a) the zero mode 
7](t) vs. t, (b) the quantum fluctuation operator ffS(r) vs. r, (c) the number of particles gN{r) 
vs. r, (d) the particle distribution gNk(r) vs. k at t = 50.8 (dashed line) and r = 399.4 (solid 
line), and (e) the ratio of the pressure and energy density p(r)/e(r) vs. r for the parameter values 
m 2 = -1, rj(T ) = 4, 7?(r ) = 0, g = 1(T 12 , h(r ) = 0.1. 

produced by parametric amplification is concentrated at higher momenta, k ~ 1. 

Figure 4 shows the corresponding case with a matter dominated background. The results 
are qualitatively very similar to those described for fig. 3 above. Due to the faster expansion, 
the zero mode (fig. 4a) finds one of the two wells more quickly and slightly less particles are 
produced. For late times, the fluctuation gT,(t) (fig.4b) decays as \/a 2 {t) oc 1/t 4 / 3 . Again we 
see an equation of state (figs. 4e) which evolves from a state between that of pure radiation 
or matter toward one of cold matter. 
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A larger Hubble constant prevents significant particle production unless the initial am- 
plitude of the zero mode is likewise increased such that the relation 77 (to) 3> M^o) is satisfied. 
For very large amplitude rj(to) ^> 1, to the extent that the mass term can be neglected and 
while the quantum fluctuation term has not grown to be large, the equations of motion 
( 4.40| ) are scale invariant with the scaling rj — > fir], H — > fiH, t —> t/ fi, and k — > [ik, where 



\i is an arbitrary scale. 

Case 3: m 2 > 0, r/(to) ^ 1- The final case we examine is that of a simple chaotic 
scenario with a positive mass term in the Lagrangian. Again, the FRW stage occurs after 
the inflationary expansion; this allows us to take zero initial temperature. 

Figure 5 shows this situation in the large N approximation for a matter dominated 
cosmology. The zero mode, f]{r), oscillates in time while decaying in amplitude from its 
initial value of r/(r = 0) = 5, r)(r = 0) = (fig.5a), while the quantum fluctuation, g£(i), 
grows rapidly for early times due to parametric resonance (figs. 5b). We choose here an 
initial condition on the zero mode which differs from that of figs 2-3 above since there is 
no significant growth of quantum fluctuations for smaller initial values. From fig.5d, we see 
that there exists a single unstable band at values of roughly k — 1 to k — 3, although careful 
examination reveals that the unstable band extends all the way to k — 0. The equation of 
state is depicted by the quantity p(r)/e(r) in fig.5e. As expected in this massive theory, the 
equation of state is matter dominated. 

First, we note that, for early times when g^(r) <C 1, the zero mode is well fit by the 
function r/(r) = r/o/(r)/a(r) where /(r) is an oscillatory function taking on values from — 1 
to 1. This is clearly seen from the envelope function r)o/a(r) shown in fig. 6a (recall that 
(?£(t) <C 1 during the entire evolution in this case). Second, the momentum that appears 



in the equations for the modes ( 4.36 ) is the physical momentum k/a(r). We therefore write 
the approximate expressions for the locations of the forbidden bands in FRW by using 
the Minkowski results of §J with the substitutions i]q — > 77/Q / a 2 (r) (where the factor of 7 
accounts for the difference in the definition of the non-linear coupling between this study 
and 0) and q 2 -> k 2 /a 2 (r). 

Making these substitutions, we find for the location in comoving momentum k of the 
forbidden band in the large N (fig. 5-6) case: 

2 

< k 2 < ^ . (6.8) 

The important feature to notice is that while the location of the unstable band (to a first 
approximation) in the case of the continuous O(N) theory is the same as in Minkowski and 
does not change in time. Again, the qualitative dynamics remains largely unchanged from 
the case of a smaller Hubble constant. 

As in the symmetry broken case of figs. 1-4, the equations of motion for large amplitude 
and relatively early times are approximately scale invariant. In fig. 6 we show the case of 
the large N evolution in a radiation dominated universe with initial Hubble constant of 
H{tq) = 2 with appropriately scaled initial value of the zero mode of t]{tq) = 16. 
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Figure 5a Figure 5b 
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FIG. 5. Symmetry unbroken, chaotic, large N, matter dominated evolution of (a) the zero 
mode ry(r) vs. r, (b) the quantum fluctuation operator gS(r) vs. r, (c) the number of particles 
gN(r) vs. r, (d) the particle distribution gNk(r) vs. k at r = 77.4 (dashed line) and r = 399.7 
(solid line), and (e) the ratio of the pressure and energy density p(r)/e(r) vs. r for the parameter 
values m 2 = +1, t?(t ) = 5, 77(70) = 0, g = 1(T 12 , H(t ) = 0.1. 



C. Late Time Behavior 

We see clearly from the numerical evolution that in the case of a symmetry broken 
potential, the late time large N solutions obey the sum rule (|5.40p . This sum rule is a 



consequence of the late time Ward identities which enforce Goldstone's Theorem. Because 
of this sum rule, we can write down the analytical expressions for the late time behavior of 
the fluctuations and the zero mode. Using eq.( |5.40D , the mode equation (|4.36| ) becomes 
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FIG. 6. Symmetry unbroken, chaotic, large N, radiation dominated evolution of (a) the zero 
mode t](t) vs. t, (b) the quantum fluctuation operator <?£(t) vs. t, (c) the number of particles 
gN(r) vs. t, (d) the particle distribution gN k {r) vs. k at r = 102.1 (dashed line) and r = 251.6 
(solid line), and (e) the ratio of the pressure and energy density p(r)/e(r) vs. r for the parameter 
values m 2 = +1, r/(r ) = 16, t?(t ) = 0, g = 10~ 12 , h{r ) = 2.0. 
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(6.9) 



This equation can be solved exactly if we assume a power law dependence for the scale factor 
a(r) = (t/to)™ with solution 



U k (r) = c k r^- 3 ")/ 2 



n-1 



+ d k r^' 2 Y,_ 



'kr^T X -'' 
n-1 



(6.10) 
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where J v and Y v are Bessel and Neumann functions respectively, and the constants c& and 
dk carry dependence on the initial conditions and on the dynamics up to the point at which 
the sum rule is satisfied. 

These functions have several important properties. In particular, in radiation or matter 
dominated universes, n < 1, and for values of wavenumber satisfying k ^> r _( - 1_n - ) /t£, the 
mode functions decay in time as l/a(r) ~ r~ n . Since the sum rule applies for late times, 
t — ro ^> 1 in dimensionless units, we see that all values of k except a very small band about 
k = redshift as l/a(r). The k = mode for any scale factor a(r) takes the form 



U {t) = C 1 + C 2 



dr' 



o?(t') 



where we used eq. ( |6.9|) and C\ and C 2 are constants depending on the initial conditions. 

We see that the k = mode freezes out for late times tending to a constant. This 
explains the support evidenced in the numerical results for values of small k (see figs. 1,3). 

These results mean that the quantum fluctuation has a late time dependence of 
(7r 2 (r)) r ~ l/a 2 (r). The late time dependence of the zero mode is given by this expres- 
sion combined with the sum rule (|5.40|) . These results are accurately reproduced by our 
numerical analysis. Note that qualitatively this late time dependence is independent of the 
choice of initial conditions for the zero mode, except that there is no growth of modes near 
k = in the case in which particles are produced via parametric amplification (figs. 4,5). 

For the radiation n = \ and matter dominated n — | universes, eq.( |6.10| ) reduces to 
elementary functions: 
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(MD) . (6.11) 



It is also of interest to examine the n > 1 case. Here, the modes of interest satisfy the 
condition k r n_1 /tQ for late times. These modes are constant in time and one sees that 
the modes are frozen. In the case of a de Sitter universe, we can formally take the limit 
n —>■ oo and we see that all modes become frozen at late times. This case is detailed in sec. 

VII p]. 



D. Discussion, Conclusions and Further Results for the FRW background 

We have shown that there can be significant particle production through quantum fluctu- 
ations after inflation [[EJ. However, this production is somewhat sensitive to the expansion 
of the universe. From our analysis of the equation of state, we see that the late time dy- 
namics is given by a matter dominated cosmology. We have also shown that the quantum 
fluctuations of the inflaton decay for late times as l/a 2 (t), while in the case of a symmetry 
broken inflationary model, the inflaton field moves to the minimum of its tree level poten- 
tial. The exception to this behavior is the case when the inflaton begins exactly at the 
unstable extremum of its potential for which the fluctuations grow out to the minimum of 
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the potential and do not decay. Initial production of particles due to parametric amplifi- 
cation is significantly greater in chaotic scenarios with symmetry broken potentials than in 
the corresponding theories with positive mass terms in the Lagrangian, given similar initial 
conditions on the zero mode of the inflaton. 



In ref. (16| we further investigate a symmetry breaking phase transition triggered by the 
lowering of the temperature as l/a(t) both in radiation dominated and matter dominated 
FRW spacetimes. We identify three different time scales: an early regime dominated by 
linear instabilities and the exponential growth of long-wavelength fluctuations, an interme- 
diate scale when the field fluctuations probe the broken symmetry states and an asymptotic 
scale wherein a scaling regime emerges for modes of wavelength comparable to or larger 
than the horizon. The scaling regime is characterized by a dynamical physical correlation 
length ^phys = dff(t) with duif) the size of the causal horizon, thus there is one correlated 
region per causal horizon. Inside these correlated regions the field fluctuations sample the 
broken symmetry states. The amplitude of the long-wavelength fluctuations becomes non- 
perturbatively large due to the early times instabilities and a semiclassical but stochastic 
description emerges in the asymptotic regime. In the scaling regime, the power spectrum is 
peaked at zero momentum revealing the onset of a Bose-Einstein condensate. The scaling 
solution results in that the equation of state of the scalar fields is the same as that of the 
background fluid. This implies a Harrison-Zeldovich spectrum of scalar density perturba- 
tions for long-wavelengths. We discuss the corrections to scaling as well as the universality 
of the scaling solution and the differences and similarities with the classical non-linear sigma 
model. 



VII. SCALAR FIELD DYNAMICS IN A FIXED INFLATIONARY 
BACKGROUND (THE DE SITTER UNIVERSE) 

We describe in this section the scalar field evolution in a fixed de Sitter background using 
the large N approximation, postponing the evolution with a dynamical background to sec. 
VII. 

We consider an initial state at a non-zero temperature Tj. This change on the initial 
conditions does not affect the initial values of the mode functions in eqs. (|5.36| ) or fl5.38|) . 

Only the expression for the quantum fluctuations E(r) changes due to the fact that the 
expectation value of the product of a creation and an annihilation operator is temperature 
dependent [[TT]. We now have instead of eq. ( |5.39| ) 
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2Tj qa(r) 



+ 
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(7.1) 



13 



We consider the case m? < with a critical temperature T c such that Tj > T c ^> |m 
The symmetry is initially unbroken and through the expansion of the universe the effective 
temperature decreases as Ti/a(t). Therefore, there is a symmetry breaking phase transition 
after a few e-folds of inflation. When the temperature falls below the critical value, the 
effective mass becomes negative. As will be seen explicitly below, when this occurs, long- 
wavelength modes become unstable and grow. Local thermodynamic equilibrium will set in 
again if the contribution from the quantum fluctuations can grow and adjust to compensate 
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for the negative mass terms on the same time scales as that in which the temperature drops. 
However, as discussed below, for very weak coupling the important time scales for the non- 
equilibrium fluctuations are of the order of [H/m 2 R ] m[l/A], which are much longer than the 
time it takes for the temperature to drop well below the critical value to practically zero. 
Thus, the non-equilibrium dynamics will proceed as if the phase transition occured via a 
quench, that is with an effective mass term, 



m 



eff 



(t) = m 2 6(ti -t)-m 2 R 6{t - U 
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± c 
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> 0. 



(7.2) 



Therefore, we choose the initial conditions on the mode functions at t% 
terms of the effective mass, 



to be given in 



Ml = m\ 



T 2 

- 1 

J>2 



A 



Ti 



(7.3) 



A. Evolution for (f)(0) = <fi(0) = 0. Analytical Results 



We begin by considering the broken symmetry situation in which the expectation value 
of the inflaton field sits atop the potential hill with zero initial velocity. This situation 
is expected to arise if the system is initially in local thermodynamic equilibrium an ini- 
tial temperature larger than the critical temperature and cools down through the critical 
temperature in the absence of an external field or bias. 

The order parameter and its time derivative vanish in the local equilibrium high tem- 
perature phase, and this condition is a fixed point of the evolution equation for the zero 
mode of the inflaton. There is no rolling of the inflaton zero mode in this case, although the 
fluctuations will grow and will be responsible for the dynamics. 

We can understand the early stages of the dynamics analytically as follows. For very 
weak coupling and early time we can neglect the backreation in the mode equations, which 
become, 



d2 ot. d 
+ 3h— + 



q 



dr 2 
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(r) 



- 1 



f 9 (r) = , 



(7.4) 



/,(0) = 



q 1 



= \jq 2 + r 2 - 1 ; r = ^ . 



(7.5) 



The solutions are of the form, 



f q (r) = exp[-^r] {a(q) J v {z) + b{q) J_ u (z)} ; z = | exp[-hr] ; v = J + ^ , (7.6) 



where the coefficients a(q) and b(q) are determined by the initial conditions: 

3 
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7r q 
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(7.7) 
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For long times, e > q/h, these mode functions grow exponentially, 



f q (r) ~ b(q) J_ u (z) 



b(q) 



2h 



T(l-u) 



(u-3/2)hr 



(7.9) 



The Bessel functions appearing in the expression for the modes f q {r) can be approxi- 
mated by their series expansion, 



f,(r) 



1/ \2 4/i 2 h 
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(7.10) 



This is an expansion in powers of q 2 /{vh 2 ) and we conclude that <?X(t) is dominated by the 
modes with q < \[h. 

The integral for g'E(r) can be approximated by keeping only the modes g < fVh, 
where / is a number of order one, and by neglecting the subtraction term which will cancel 
the contributions from high momenta. Numerically, even with the backreaction taken into 
account, the integral is dominated by modes q < f ~ 10 — 20 in all of the cases that we 



studied (see ref. [13 



The contribution to the fluctuations from these unstable modes is: 



J f 3 fr 3/2 rm\ ( Ml 
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(7.11) 



where again, we have taken the high temperature limit, Tj ~ T c ^> rriR. 

From this equation, we can estimate the value of r s , the 'spinodal time', at which the 
contribution of the quantum fluctuations becomes comparable to the contribution from the 
tree level terms in the equations of motion. This time scale is obtained from the condition 
gZ(r s ) = 0(1): 



1 



(2u - 3)h 
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(7.12) 



which is in good agreement with our numerical results (see ref. [0)- F° r values of h > 1, 
which, as argued below, lead to the most interesting case, an estimate for the spinodal time 
is, 



3h 



ln[l/^]+C(l) 



(7.13) 



which is consistent with our numerical results (see [|i~3[| ). 

For r > t s , the effects of backreaction become very important, and the contribution from 
the quantum fluctuations competes with the tree level terms in the equations of motion, 
shutting-off the instabilities. Beyond r s , only a full numerical analysis will capture the 
correct dynamics. 

It is worth mentioning that had we chosen zero temperature initial conditions, then the 
coupling g — > g (see fig.0) and the estimate for the spinodal time would have been, 
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r s ~— ln[l/0] + 0(l), 



(7.14) 



that is, roughly a factor 2 larger than the estimate for which the de Sitter stage began at 
a temperature above the critical value. Therefore, eq. ( [7.130 represents an underestimate of 
the spinodal time scale at which fluctuations become comparable to tree level contributions. 

The number of e-folds occurring during the stage of growth of spinodal fluctuations is 
therefore, 

3h 2 3h 2 
M e ~ —^-\n[l/y/g\ (high temperature) and Af e ~ —— hi[l/<?] (zero temperature) 

It is a factor 2 larger for zero temperature. Thus, it becomes clear that with g m 1CT 12 
and h > 2, a required number of e-folds, M e ~ 100 can easily be accommodated before the 
fluctuations become large, modifying the dynamics and the equation of state. 

The implications of these estimates are important. The first conclusion drawn from these 
estimates is that a 'quench' approximation is well justified (see [|13|)- While the temperature 
drops from an initial value of a few times the critical temperature to below critical in just 
a few e-folds, the contribution of the quantum fluctuations needs a large number of e-folds 
to grow to compensate for the tree-level terms and overcome the instabilities. Only for a 
strongly coupled theory is the time scale for the quantum fluctuations to grow short enough 
to restore local thermodynamic equilibrium during the transition. 

The second conclusion is that most of the growth of spinodal fluctuations occurs during 
the inflationary stage, and with g « 10~ 12 and H > thr, the quantum fluctuations become 
of the order of the tree-level contributions to the equations of motion within the number of 
e-folds necessary to solve the horizon and flatness problems. Since the fluctuations grow to 
become of the order of the tree level contributions at times of the order of this time scale, 
for larger times they will modify the equation of state substantially and will be shown in 
sec. IX to provide a graceful exit from the inflationary phase within an acceptable number 
of e-folds. 



B. The late time limit 



For late the dynamics freezes out. The fluctuation, ^S(t) = 1, and the mode functions 
effectively describe free, minimally coupled, massless particles. The sum rule, 



1 + #£(oo) = 0, 



(7.15) 



is obeyed exactly in the large N limit as in the Minkowski case [BLPfl. 

We now show that this value is a self-consistent solution of the equations of motion for 
the mode functions, and the only stationary solution for asymptotically long times. 

In the late time limit, the effective time dependent mass term, — 1 + rf + gH, in the 



equation for the mode functions, ( 5.371) , vanishes (in this case with r\ = 0). Therefore, these 
mode equations asymptotically become, 
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The general solutions are given by, 
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where Jspi(z) an d ^3/2 (-2) are the Bessel and Neumann functions, respectively. The coeffi- 
cients, c±(q) can be computed for large q by matching fq Sy (r) with the WKB approximation 
to the exact mode functions f q (r) that obey the initial conditions ( 5.38p . The WKB ap- 

and we find for large q, 



proximation to f q (r) has been computed in ref. [p 



c±(q) 
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(7.19) 



In the t — > 00 limit, we have for fixed q, 




c-{q) ■ 



(7.20) 



which are independent of time asymptotically, and explains why the power spectrum of 
quantum fluctuations freezes at times larger than the spinodal. This behavior is confirmed 



numerically [see ref. [13]]. Clearly at early times the mode functions grow exponentially, 
and at times of the order of r s , when <?£(t) w 1 the mode functions freeze-out and become 
independent of time. Notice that the largest q modes have grown the least, explaining why 
the integral is dominated by q < 10 — 20. 

For asymptotically large times, ^S(r) is given by, 



pE(oo) = g h 2 [ 
Jo 



00 * coth ( ^ 
q \2T 
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(7.21) 



where only one term in the UV subtraction survived in the r = 00 limit. The factor coth f |jH 
in eq. ([7.21[) takes into account the nonzero initial temperature T. 

For consistency, this integral must converge and be equal to 1 as given by the sum rule. 
For this to be the case and to avoid the potential infrared divergence in ( |7.21|) , the coefficients 
c_(g) must vanish at q = 0. The mode functions are finite in the q — > limit provided, 



c_(g) 9 =° C q 



3/2 



(7.22) 



where C is a constant. 

The numerical analysis clearly shows that the mode functions remain finite as q — > 0, 
and the coefficient C can be read off from these figures. This is a remarkable result. It is well 
known that for free massless minimally coupled fields in de Sitter space-time with Bunch- 
Davies boundary conditions, the fluctuation contribution (n 2 (x,t)) grows linearly in time 
as a consequence of the logarithmic divergence in the integrals ||28|| . However, in our case, 
although the asymptotic mode functions are free, the coefficients that multiply the Bessel 
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functions of order 3/2 have all the information of the interaction and initial conditions 
and must lead to the consistency of the sum rule. Clearly the sum rule and the initial 
conditions for the mode functions prevent the coefficients c±(q) from describing the Bunch- 
Davies vacuum. These coefficients are completely determined by the initial conditions and 
the dynamics. This is the reason why the fluctuation freezes at long times unlike in the free 
case in which they grow linearly [p8fl . 

It is easy to see from eqs.( [5.42"D and ( |7.20| ) that the energy and pressure vanish for r — > oo. 



C. Discussion and Conclusions for the de Sitter background 



We have identified analytically and numerically two distinct regimes for the dynamics 
determined by the initial condition on the expectation value of the zero mode of the inflaton 

HI- 

1. When i](0) « g 1 / 4 (or g 1 ! 2 for Tj = 0), the dynamics is driven by quantum (and 
thermal) fluctuations. Spinodal instabilities grow and eventually compete with tree 
level terms at a time scale, r s > — 3h\n[g}/2. The growth of spinodal fluctuations 
translates into the growth of spatially correlated domains which attain a maximum 
correlation length (domain size) of the order of the horizon. For very weak coupling 
and h > 1 this time scale can easily accommodate enough e-folds for inflation to solve 
the flatness and horizon problems. The quantum fluctuations modify the equation 
of state dramatically providing a means for a graceful exit to the inflationary stage 
without slow-roll. 

This non-perturbative description of the non-equilibrium effects in this regime in which 
quantum (and thermal) fluctuations are most important provides a reliable under- 
standing of the relevant non-perturbative, non-equilibrium effects of the fluctuations 
that have not been revealed before in this setting || - |]I4| . 

These initial conditions are rather natural if the de Sitter era arises during a phase 
transition from a radiation dominated high temperature phase in local thermodynamic 
equilibrium, in which the order parameter and its time derivative vanish. 

2. When r](0) >> g l l A (or g x l 2 for Tj = 0), the dynamics is driven solely by the classical 
evolution of the inflaton zero mode. The quantum and thermal fluctuations are always 
perturbatively small (after renormalization) , and their contribution to the dynamics is 
negligible for weak couplings. The de Sitter era will end when the kinetic contribution 
to the energy becomes of the same order as the 'vacuum' term. This is the realm 
of the slow-roll analysis whose characteristics and consequences have been analyzed 
in the literature at length. These initial conditions, however, necessarily imply some 
initial state either with a biasing field that favors a non-zero initial expectation value, 
or that in the radiation dominated stage, prior to the phase transition, the state was 
strongly out of equilibrium with an expectation value of the zero mode different from 
zero. Although such a state cannot be ruled out and would naturally arise in chaotic 
scenarios, the description of the phase transition in this case requires further input on 
the nature of the state prior to the phase transition. 
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VIII. SELF-CONSISTENT EVOLUTION OF MATTER FIELDS WITH A 
DYNAMICAL COSMOLOGICAL BACKGROUND 



We present in this section the full self-consistent matter-geometry dynamics |[L4|| . That 



is, the scale factor a{t) is here a dynamical variable determined by the Einstein-Friedman 
eq. (|5.41|) - (|5.42|) coupled with the scalar field evolution eqs.( |5.37D - (|5.38|) . 

In order to provide the full solution we now must provide the values of i](0), r](0), and 
ho. Assuming that the inflationary epoch is associated with a phase transition at the GUT 
scale, this requires that Nm A R / g pa (10 15 Gev ) 4 and assuming the bound on the scalar self- 
coupling g pa 10~ 12 — 10~ 14 (this will be seen later to be a compatible requirement), we find 
that ho ~ iV 1 / 4 which we will take to be reasonably given by h pa 1 — 10 (for example in 
popular GUT's iV pa 20 depending on particular representations). 

We will begin by studying the case of most interest from the point of view of describing 
the phase transition: 7/(0) = and 7/(0) = 0, which are the initial conditions that led to 
puzzling questions. With these initial conditions, the evolution equation for the zero mode 
eq. (|5.371) determines that T)(t) = by symmetry. 



A. Early time dynamics: 

Before engaging in the numerical study, it proves illuminating to obtain an estimate of 
the relevant time scales and an intuitive idea of the main features of the dynamics. Because 
the coupling is so weak (g ~ 10~ 12 <C 1) and after renormalization the contribution from 
the quantum fluctuations to the equations of motion is finite, we can neglect all the terms 
proportional to g in eqs. (|5.42[) and ( p.37| ). 



For the case where we choose T](t) = and the evolution equations for the mode functions 
are those for an inverted oscillator in De Sitter space-time, which have been studied in sec. 



VII [[H]]. One obtains the approximate solutions (|7.6|) -(|7T 

After the physical wavevectors cross the horizon, i.e. when qe~ haT /ho < 1 we find that 
the mode functions factorize: 



m qe ~ h = <<hQ f^nj e (»-m h0 r. (8 . 1} 

This result reveals a very important feature: because of the negative mass squared term 
in the matter Lagrangian leading to symmetry breaking (and v > 3/2), we see that all 
of the mode functions grow exponentially after horizon crossing (for positive mass squared 
v < 3/2, and they would decrease exponentially after horizon crossing). This exponential 
growth is a consequence of the spinodal instabilities which is a hallmark of the process of 
phase separation that occurs to complete the phase transition. We note, in addition that 
the time dependence is exactly given by that of the q = mode, i.e. the zero mode, which is 
a consequence of the redshifting of the wavevectors and the fact that after horizon crossing 
the contribution of the term q 2 /a 2 (r) in the equations of motion become negligible. We 
clearly see that the quantum fluctuations grow exponentially and they will begin to be of 



V 
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the order of the tree level terms in the equations of motion when gTt(r) 
substantially before the end of inflation 



1. At large times 



K e 



(2v-3)h T 



with J-'(ho) a finite constant that depends on the initial conditions and is found numerically 
to be of 0(1) [see fig.0. 

In terms of the initial dimensionful variables, the condition gH{r) ~ 1 translates to 
< 7i 2 (x, t) >/?~ 2m 2 R /g, i.e. the quantum fluctuations sample the minima of the (renormal- 
ized) tree level potential. We find that the time at which the contribution of the quantum 
fluctuations becomes of the same order as the tree level terms is estimated to be [|13[ 
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(8.2) 



At this time, the contribution of the quantum fluctuations makes the back reaction very im- 
portant and, as will be seen numerically, this translates into the fact that r s also determines 
the end of the De Sitter era and the end of inflation. The total number of e-folds during the 
stage of exponential expansion of the scale factor (constant ho) is given by 
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(8.3) 



For large ho we see that the number of e-folds scales as h\ as well as with the logarithm of 



the inverse coupling. These results ( |8.1| - |8.3| ) will be confirmed numerically below and will be 
of paramount importance for the interpretation of the main consequences of the dynamical 
evolution. 

As discussed in sec. VII. C, the early time dynamics is dominated by classical or quantum 
effects depending on the ratio between the time scales r c and r s . 

If r c is much smaller than the spinodal time t s given by eq.( |8.2|) then the classical 
evolution of the zero mode will dominate the dynamics and the quantum fluctuations will 
not become very large, although they will still undergo spinodal growth. On the other hand, 
if r c ^> t s the quantum fluctuations will grow to be very large well before the zero mode 
reaches the non-linear regime. In this case the dynamics will be determined completely by 
the quantum fluctuations. Then the criterion for the classical or quantum dynamics is given 
by 



r]{0) ^>V9~h 



classical dynamics 
quantum dynamics 



(8.4) 



or in terms of dimensionful variables 0(0) ^> H leads to classical dynamics and 0(0) -C H 
leads to quantum dynamics. 

However, even when the classical evolution of the zero mode dominates the dynamics, 
the quantum fluctuations grow exponentially after horizon crossing unless the value of <fi(t) 
is very close to the minimum of the tree level potential. In the large N approximation the 
spinodal line, that is the values of <p(t) for which there are spinodal instabilities, reaches all 
the way to the minimum of the tree level potential as can be seen from the equations of 
motion for the mode functions. Therefore even in the classical case one must understand 
how to deal with quantum fluctuations that grow after horizon crossing. 
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B. Numerics 



The time evolution is carried out by means of a fourth order Runge-Kutta routine with 
adaptive step-sizing while the momentum integrals are carried out using an 11-point Newton- 
Cotes integrator. The relative errors in both the differential equation and the integration 
are of order 10~ 8 . We find that the energy is covariantly conserved throughout the evolution 
to better than a part in a thousand. Figs. show gT,(r) vs. r, h{r) vs. r and In \f q {j)\ 2 
vs. t for several values of q with larger q's corresponding to successively lower curves. Figs. 
10 show p(t)/s(t) and the horizon size /i _1 (r) for g = 10 
we have chosen the representative value ho = 2.0. 
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FIG. 7. flE(r) vs. r, for r?(0) = 0,r)(0) = 0, g = 10~ 14 , h = 2.0 



FIG. 8. h(r) vs. r, for 77(0) = 0, r)(0) = 0, 5 = 10" 14 , /i = 2.0. 

Figs. and § show clearly that when the contribution of the quantum fluctuations <?£(t) 
becomes of order 1 inflation ends, and the time scale for S'S(t) to reach 0(1) is very well 
described by the estimate (|8.2|) . From fig.[7] we see that this happens for r = t s w 90, leading 
to a number of e-folds N e m 180 which is correctly estimated by eqs. ( |8.2j )-(|8~^). 

Fig. |^ shows clearly the factorization of the modes after they cross the horizon as 
described by eq.( |8.1| ). The slopes of all the curves after they become straight lines in fig.|9] is 
given exactly by (2u— 3), whereas the intercept depends on the initial condition on the mode 
function and the larger the value of q the smaller the intercept because the amplitude of the 
mode function is smaller initially. Although the intercept depends on the initial conditions 
on the long-wavelength modes, the slope is independent of the value of q and is the same as 
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what would be obtained in the linear approximation for the square of the zero mode at times 
long enough that the decaying solution can be neglected but short enough that the effect 
of the non-linearities is very small. Notice from the figure that when inflation ends and the 
non-linearities become important all of the modes effectively saturate. This is also what one 
would expect from the solution of the zero mode: exponential growth in early-intermediate 
times (neglecting the decaying solution), with a growth exponent given by {y — 3/2) and an 
asymptotic behavior of small oscillations around the equilibrium position, which for the zero 
mode is 77 = 1, but for the q 7^ modes depends on the initial conditions. All of the mode 
functions have this behavior once they cross the horizon. We have also studied the phases 
of the mode functions and we found that they freeze after horizon crossing in the sense that 
they become independent of time. This is natural since both the real and imaginary parts 
of f q (r) obey the same equation but with different boundary conditions. After the physical 
wavelength crosses the horizon, the dynamics is insensitive to the value of q for real and 
imaginary parts and the phases become independent of time. Again, this is a consequence 
of the factorization of the modes. 




FIG. 9. ln|/ g (r)| 2 vs. r, for 7/(0) = 0,j)(0) = 0,g = 10 
with smaller q corresponding to larger values of In |/q(r)| 2 . 



-14 



,ho = 2.0 for q = 0.0,5,10,15,20 



The growth of the quantum fluctuations is sufficient to end inflation at a time given by 
t s in eq. (|8.2| ). Furthermore fig. [H] shows that during the inflationary epoch p(r)/e{r) ~ — 1 
and the end of inflation is rather sharp at r s with p(r)/e(r) oscillating between ±1 with 



zero average over the cycles, resulting in matter domination. Fig. [11] shows this feature very 
clearly; h(r) is constant during the de Sitter epoch and becomes matter dominated after 



the end of inflation with h (r) 



There are small oscillations around this value 



because both p(r) and e(r) oscillate. These oscillations are a result of small oscillations 
of the mode functions after they saturate, and are also a feature of the solution for a zero 
mode. 

All of these features hold for a variety of initial conditions. As an example, we show in 
ref. H the case of an initial Hubble parameter of ho = 10. 
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FIG. 10. p/e vs. r, for 77(0) = 0, t)(0) = 0, g = 10~ 14 , h = 2.0. 
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FIG. 11. l/h(r) vs. r, for r/(0) = 0, j)(0) = 0,5 = 10~ 14 , /i = 2.0. 

C. Zero Mode Assembly: 

This remarkable feature of factorization of the mode functions after horizon crossing can 
be elegantly summarized as 

fk(t)\ kvh(t )^H = g{q, h)fo(r), (8.5) 

with k p h(t) = ke~ Ht being the physical momentum, g(q, h) a complex constant, and /o( r ) a 
real function of time that satisfies the mode equation with q = and real initial conditions 
which will be inferred later. 

Then we consider the contribution of these modes to the renormalized quantum fluctu- 
ations a long time after the beginning of inflation (so as to neglect the decaying solutions), 
we find that 

#£(r) w Ce {2u - 3)hT + small , 

where 'small' stands for the contribution of mode functions associated with momenta that 
have not yet crossed the horizon at time r, which give a perturbatively small (of order g) 
contribution. We find that several e-folds after the beginning of inflation and till inflation 
ends, this factorization of superhorizon modes implies the following: 

gJq 2 dq\f*(T)\^\C ff (T) , 9 Jq 2 dq\fi(T)\^\C ffZ(T), 
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9 



a 2 (r 



-dq \f (r) 



a 2 (r 



(8.6) 



where we have neglected the weak time dependence arising from the perturbatively small 
contributions of the short-wavelength modes that have not yet crossed the horizon, and the 
integrals above are to be understood as the fully renormalized (subtracted), finite integrals. 
For rj = 0, we note that (|8.6j) and the fact that /o(t) obeys the equation of motion for the 
mode with q = leads at once to the conclusion that in this regime [<?£(t)] 5 = \C \f (T) 
obeys the zero mode equation of motion 



^ + 3^-l + (|C |/o(r)) 2 



|C |/o(r) = 



(8.7) 



It is clear that at a time T4 several e-folds after the beginning of inflation, we can define an 
effective zero mode as 



^ 2 //( r ) 



<?E(t), or in dimensionful variables, 



(t) 



7T 



(£,*)) 



R 



(8.8) 



Although this identification seems natural, we emphasize that it is by no means a trivial or 
ad- hoc statement. There are several important features that allow an unambiguous identifi- 
cation: i) [{ir 2 (x, t))u] is a fully renormalized operator product and hence finite, ii) because 
of the factorization of the superhorizon modes that enter in the evaluation of [(tt 2 (x, t))u], 
(f) e ff(t) (|8.8j ) obeys the equation of motion for the zero mode, iii) this identification is valid 
several e-folds after the beginning of inflation, after the transient decaying solutions have 
died away and the integral in (7r 2 (x,t)) is dominated by the modes with wavevector k that 
have crossed the horizon at t(k) <C t. Numerically we see that this identification holds 
throughout the dynamics except for a very few e-folds at the beginning of inflation. This 
factorization determines at once the initial conditions of the effective zero mode that can be 
extracted numerically: after the first few e-folds and long before the end of inflation we find 



<p eff (t) = <f) eff (t 



A) e 



iy- 



•|)H(t-tA) 



(8.9) 



where we parameterized 



W/(*a) = |- HH/m) 

Z7T 



to make contact with the literature. As is shown in fig. (|T^D, we find numerically that 
TiHjm) w 0(1) for a large range of 0.1 < H/m < 50 and that this quantity depends on 
the initial conditions of the long wavelength modes. 

Therefore, in summary, the effective composite zero mode obeys 



dr 2 



Veff(r) = ; Veff(r A ) 



[8.10) 



where ta is the time at which the composite mode becomes effective, Veff{ T A) = 

4> e ff(tA) is obtained numerically for a given ho by fitting the intermediate time be- 
havior of gTi(r) with the growing zero mode solution. Recall that \r = 871" 2 g. 
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FIG. 12. F(H/m) vs. h, where F(H/m) is defined by the relation 4> e ff{t A ) = (H /2ir)F(H jra) 
(see eqs. (O) and (IO)). 
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FIG. 13. r)l n (r) vs. r, for %// (0) = 3.94 x 1(T 7 , r) e// (0) = 0.317^/(0), g = 10" 14 , h 
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The initial conditions were obtained by fitting the intermediate time regime of gS(r) in fig.^. 



Veff(i~) is the solution of eq.( 8.10| ) with these initial conditions. 
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FIG. 14. h(r) vs. r, obtained from the solution of eqs. ( 8.10j ) and ( |S.12| ) with the conditions 



of fig. 13 



Furthermore, this analysis shows that in the case r] = 0, the renormalized energy and 
pressure in this regime in which the renormalized integrals are dominated by the superhori- 
zon modes are given by 



2Nm%fl + l/_ 1+r?2 ^ 



A, 



2Nm 



A, 



(8.11) 
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where we have neglected the contribution proportional to l/a 2 (r) because it is effectively 
red-shifted away after just a few e-folds. We found numerically that this term is negligible 
after the interval of time necessary for the superhorizon modes to dominate the contribution 
to the integrals. Then the dynamics of the scale factor is given by 



h 2 



+ 4(- 1 +^// 



f8.12) 



We have numerically evolved the set of effective equations (8.1U)-( 8.12|) by extracting the 
initial condition for the effective zero mode from the intermediate time behavior of g£(r). 
We found a remarkable agreement between the evolution of r\ 2 e ^ and g^(r) and between the 
dynamics of the scale factor in terms of the evolution of ^//(t), and the full dynamics of 
the scale factor and quantum fluctuations within our numerical accuracy. Figs. [13| and |14] 
show the evolution of ^e//( T ) an d h(r) respectively from the classical evolution eqs. ( |3.10| ) 
and 12 ) using the initial condition ^//(t^) extracted from the exponential fit of gSfV) in 
the intermediate regime. These figures should be compared to figs. [7| and |8|. We have also 
numerically compared p/e given solely by the dynamics of the effective zero mode and it is 
again numerically indistinguishable from that obtained with the full evolution of the mode 
functions. 

This is one of the main results of our work 



14]. In summary: the modes that become 



superhorizon sized and grow through the spinodal instabilities assemble themselves into an 
effective composite zero mode a few e-folds after the beginning of inflation. This effective 
zero mode drives the dynamics of the FRW scale factor, terminating inflation when the non- 
linearities become important. In terms of the underlying fluctuations, the spinodal growth of 
superhorizon modes gives a non-perturbatively large contribution to the energy momentum 
tensor that drives the dynamics of the scale factor. Inflation terminates when the mean 
square root fluctuation probes the equilibrium minima of the tree level potential. 

The extension of this analysis to the case for which r](0) ^ is straightforward. Since 
both 7](t) and JgT^r) = |C |/o( r ) obey the equation for the zero mode, eq.( |5.37| ), it is clear 
that we can generalize our definition of the effective zero mode to be 



Veff (r) = ^ 2 (r)+pE(r) . 
which obeys the equation of motion of a classical zero mode: 



(8.13) 



^ + 3/ i (r)A_ 1 + , e//(r) 



Ve ff (r) = . (8.14) 
If this effective zero mode is to drive the FRW expansion, then the additional condition 



Veff fo - 2 Veff Veff fo fo + vlff fo = Veff fo - Veff fo 







(8.15) 



must also be satisfied. One can easily show that this relation is indeed satisfied if the 
mode functions factorize as in ( |3.5| ) and if the integrals in eqs. (|8.6| ) are dominated by the 
contributions of the superhorizon mode functions. This leads to the conclusion that the 
gravitational dynamics is given by eqs. ( |8.11| ) - ( |8.12| ) with r] e ff(r) defined by ( ^.13| ). 
Eq.( |S.15| ) is just the vanishing of the wronskian of r] e ff(r) and /o( r )- Namely, rj e ff(r) and 
/ (t) just differ in a constant factor. 
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We see that in all cases, the full large N quantum dynamics in these models of inflation- 
ary phase transitions is well approximated by the equivalent dynamics of a homogeneous, 
classical scalar field with initial conditions on the effective field ^//(ta) — \[9 h J-'(ho). 
We have verified these results numerically for the field and scale factor dynamics, finding 
that the effective classical dynamics reproduces the results of the full dynamics to within 
our numerical accuracy. We have also checked numerically that the estimate for the classical 
to quantum crossover given by eq.( fj.4[ ) is quantitatively correct. Thus in the classical case 
in which r/(0) ^> ^fg h we find that r] e ff(r) = r](r), whereas in the opposite, quantum case 

Veff(r) = yJgZ(r). 

This remarkable feature of zero mode assembly of long-wavelength, spinodally unstable 
modes is a consequence of the presence of the horizon. It also explains why, despite the fact 
that asymptotically the fluctuations sample the broken symmetry state, the equation of state 
is that of matter. Since the excitations in the broken symmetry state are massless Goldstone 
bosons one would expect radiation domination. However, the assembly phenomenon, i.e. 
the redshifting of the wave vectors, makes these modes behave exactly like zero momentum 
modes that give an equation of state of matter (upon averaging over the small oscillations 
around the minimum). 

Subhorizon modes at the end of inflation with q > ho e h ° Ts do not participate in the 
zero mode assembly. The behavior of such modes do depend on q after the end of inflation. 
Notice that these modes have extremely large comoving q since h e hoTs > 10 26 . As discussed 
in sec. VI such modes decrease with time after inflation as ~ V a ( r ) EH- 

D. Making sense of small fluctuations: 

Having recognized the effective classical variable that can be interpreted as the compo- 
nent of the field that drives the FRW background and rolls down the classical potential hill, 
we want to recognize unambiguously the small fluctuations. We have argued above that 
after horizon crossing, all of the mode functions evolve proportionally to the zero mode, and 
the question arises: which modes are assembled into the effective zero mode whose dynamics 
drives the evolution of the FRW scale factor and which modes are treated as perturbations? 
In principle every k ^ mode provides some spatial inhomogeneity, and assembling these 
into an effective homogeneous zero mode seems in principle to do away with the very in- 
homogeneities that one wants to study. However, scales of cosmological importance today 
first crossed the horizon during the last 60 or so e-folds of inflation. Recently Grishchuk 
has argued that the sensitivity of the measurements of AT/T probe inhomogeneities 
on scales ~ 500 times the size of the present horizon. Therefore scales that are larger than 
these and that have first crossed the horizon much earlier than the last 60 e-folds of inflation 
are unobservable today and can be treated as an effective homogeneous component, whereas 
the scales that can be probed experimentally via the CMB inhomogeneities today must be 
treated separately as part of the inhomogeneous perturbations of the CMB. 

Thus a consistent description of the dynamics in terms of an effective zero mode plus 
'small' quantum fluctuations can be given provided the following requirements are met: 
a) the total number of e-folds N e 3> 60, b) all the modes that have crossed the horizon 
before the last 60-65 e-folds are assembled into an effective classical zero mode via e //(t) = 
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[0o(t) + (7r 2 (x, t))^] 1 , c) the modes that cross the horizon during the last 60-65 e-folds 
are accounted as 'small' perturbations. The reason for the requirement a) is that in the 
separation (p(x,t) = (fi e ff(t) + 6<p(x,t) one requires that 5(fi(x,t)/(j) e ff(t) <C 1. As argued 
above, after the modes cross the horizon, the ratio of amplitudes of the mode functions 
remains constant and given by e^ - ^^ with AN being the number of e-folds between the 
crossing of the smaller k and the crossing of the larger k. Then for 5<p(x, t) to be much smaller 
than the effective zero mode, it must be that the Fourier components of 5<f) correspond to 
very large fc's at the beginning of inflation, so that the effective zero mode can grow for a 
long time before the components of 50 begin to grow under the spinodal instabilities. In 
fact requirement a) is not very severe; in the figs.0-|IT]we have taken ho = 2.0 which is a very 
moderate value and yet for g = 10~ 12 the inflationary stage lasts for well over 100 e-folds, 
and as argued above, the larger ho for fixed g, the longer is the inflationary stage. Therefore 
under this set of conditions, the classical dynamics of the effective zero mode <fr e ff(t) drives 
the FRW background, whereas the inhomogeneous fluctuations S(f)(x,t), which are made 
up of Fourier components with wavelengths that are much smaller than the horizon at the 
beginning of inflation and that cross the horizon during the last 60 e-folds, provide the 
inhomogeneities that seed density perturbations. 



E. Scalar Metric Perturbations: 



Having identified the effective zero mode and the 'small perturbations', we are now in 
position to provide an estimate for the amplitude and spectrum of scalar metric perturba- 
tions. We use the clear formulation in ref. |37| in terms of gauge invariant variables. In 
particular we focus on the dynamics of the Bardeen potential [33], which in longitudinal 



gauge is identified with the Newtonian potential. The equation of motion for the Fourier 
components (in terms of comoving wavevectors) for this variable in terms of the effective 
zero mode is [eq.(4.48) in [[HJ 



H(t) 



4>eff(t) 
<Peff(t). 



a W \ <Peff(t). 



= o. 



(8.16) 



We are interested in determining the dynamics of for those wavevectors that cross the 
horizon during the last 60 e-folds before the end of inflation. During the inflationary stage 
and substantially before its end the numerical analysis yields to a very good approximation 



H(t) = Ho ; e// (t) = e// (t A ) e (-D ff o(*-U) ) 
where Ho is the value of the Hubble constant during inflation, leading to 



^v-2)H (t-t A ) 



'k e -H (t-t A y 



+ b k H^ 



' f~ e -H {t-t A )- 
. ^0 



(8.17) 



(8.18) 



The coefficients a^, bk are determined by the initial conditions and v is given by eq.( [7.6|) . 

Since we are interested in the wavevectors that cross the horizon during the last 60 
e-folds, the consistency for the zero mode assembly and the interpretation of 'small pertur- 
bations' requires that there must be many e-folds before the last 60. We are then considering 
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wavevectors that were deep inside the horizon at the onset of inflation. $&(£) is related to 
the canonical 'velocity field' that determines scalar perturbations of the metric and which is 
quantized with Bunch-Davies initial conditions for the large fc-mode functions. The relation 
between and v and the initial conditions on v lead at once to a determination of the 



coefficients and for k >> H [see eqs.(13.5) and (13.9) in ref. |J7 



Thus we find that the amplitude of scalar metric perturbations after horizon crossing is 
given by 



• 3 

= ki\*„{t)\ =^r(.-l) (^tP e (^^(*-*x) . (8 .20) 

7T Mp| \ k J 



The power spectrum per logarithmic k interval is given by |<5fc(t)| 2 . The time dependence 
of |<$fc(i)| displays the unstable growth associated with the spinodal instabilities of super- 
horizon modes and is a hallmark of the phase transition. This time dependence can be also 
understood from the constraint equation that relates the Bardeen potential to the gauge 
invariant field fluctuations [[37]], which in longitudinal gauge are identified with 5<p(x,t). 
The constraint equation and the evolution equations for the gauge invariant scalar field 
fluctuations are JJ7[ 



d 4:71 ■ • 

= uf a W &// . ( 8 - 21 ) 



dV ' M' Pl 



— + 3H— + — +M 2 



dt 2 



dt 



+ 2V'( 



>eff. 



(8.22) 



Since the right hand side of (|8.21|) is proportional to (f) e ff/Mp l -C 1 during the inflationary 
epoch in this model, we can neglect the terms proportional to and on the left hand 
side of (|8.22|) , in which case the equation for the gauge invariant scalar field fluctuation is 
the same as for the mode functions. In fact, since 8<jf£ is gauge invariant we can evaluate 
it in the longitudinal gauge wherein it is identified with the mode functions fk{t). Then 
absorbing a constant of integration in the initial conditions for the Bardeen variable, we find 



4tt 
M 2 Pl a(t) 



'to 



0(0 <f>eff{t') fk(t') dt' + O Uj^ 



(8.23) 



and using that (f>(t) oc e ( u - 3 / 2 ) H ot anc j that after horizon crossing fk(t) cx e ( u - 3 / 2 ) H ot ^ one 
obtains at once the time dependence of the Bardeen variable after horizon crossing. In 
particular the time dependence is found to be oc e ( 2u - 3 ) H o t _ it i s then clear that the time 
dependence is a reflection of the spinodal (unstable) growth of the superhorizon field fluc- 
tuations. 

To obtain the amplitude and spectrum of density perturbations at second horizon crossing 



we use the conservation law associated with the gauge invariant variable [37| 



53 



which is valid after horizon crossing of the mode with wavevector k. Although this conserva- 



tion law is an exact statement of superhorizon mode solutions of eq. fl8.16|) , we have obtained 



solutions assuming that during the inflationary stage H is constant and have neglected the 
H term in Eq. ( 8.16f) . Since during the inflationary stage, 



^ U \ . ZT2 (dVeffK M u-1 

l Pl 



H(t) = - w ti ff {t)<xHl < <i? o (8-25) 



and (f)/cj) ~ Ho, the above approximation is justified. We then see that (fy^fJt) oc e^ -3 ^ * 
which is the same time dependence as that of Thus the term proportional to 1/(1 + 

p/e) in Eq. ( J8.24 ) is indeed constant in time after horizon crossing. On the other hand, 



the term that does not have this denominator evolves in time but is of order (1 + p/e) 
—2H/3H 2 <C 1 with respect to the constant term and therefore can be neglected. Thus, we 
confirm that the variable £ is conserved up to the small term proportional to (1 + p/e)§k 
which is negligible during the inflationary stage. This small time dependence is consistent 
with the fact that we neglected the H term in the equation of motion for The validity 

of the conservation law has been recently studied and confirmed in different contexts ||39|| . 
Notice that we do not have to assume that vanishes, which in fact does not occur. 

However, upon second horizon crossing it is straightforward to see that $fc(t/) ~ 0. The 
reason for this assertion can be seen as follows: eq.(^.22|) shows that at long times, when 



the effective zero mode is oscillating around the minimum of the potential with a very small 
amplitude and when the time dependence of the fluctuations has saturated (see fig|J), 
will redshift as ~ l/a(t) Jl2] and its derivative becomes extremely small. 

Using this conservation law, assuming matter domination at second horizon crossing, 
and « |J7|], we find 



3 I» (2H 

Ok{tf)=—- 3 —— , 8.26 

5 7r [v - |J ^[Ho/m) V k ' 



where ^(Hq/tti) determines the initial amplitude of the effective zero mode ( |8.9|) . We can 
now read the power spectrum per logarithmic k interval 

V 8 {k) = \5 k \ 2 oc k-^~l\ (8.27) 

leading to the index for scalar density perturbations 

n s = l-2(v-^j . (8.28) 

For Ho/m ^> 1, we can expand v — 3/2 as a series in m 2 /Hq in eq.( |3.26| ). Given that 
the comoving wavenumber of the mode which crosses the horizon n e-folds before the end of 
inflation is k = H e^ Ne ~ n ^ where N e is given by ( |8.3| ), we arrive at the following expression 
for the amplitude of fluctuations on the scale corresponding to n in terms of the De Sitter 
Hubble constant and the coupling A = 87T 2 g: 
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Mt f )\ 



9H 3 



10(2vr) 3 / 2 m 3 



2m 2 



3# 2 V6 



7 7 
-ln2-^ 



m 



(8.29) 



Here, 7 is Euler's constant. Note the explicit dependence of the amplitude of density per- 
turbations on y/g. For n ~ 60, the factor exp(nm 2 /SHq) is 0(100) for H /m = 2, while it 
is 0(1) for H /m > 4. Notice that for H /m large, the amplitude increases approximately 
as (Ho/m) 3 , which will place strong restrictions on g in such models. 

We remark that we have not included the small corrections to the dynamics of the 
effective zero mode and the scale factor arising from the non-linearities. We have found 
numerically that these nonlinearities are only significant for the modes that cross about 60 
e-folds before the end of inflation for values of the Hubble parameter H /mji > 5. The effect 
of these non-linearities in the large N limit is to slow somewhat the exponential growth of 
these modes, with the result of shifting the power spectrum closer to an exact Harrison- 
Zeldovich spectrum with n s = 1. Since for Hq/itlh > 5 the power spectrum given by ( |3.28| ) 
differs from one by at most a few percent, the effects of the non-linearities are expected to 
be observationally unimportant. The spectrum given by ( 8.26| ) is similar to that obtained 
in references [j3|[3ll although the amplitude differs from that obtained there. In addition, we 
do not assume slow roll for which {y — |) <C 1, although this would be the case if N e 3> 60. 

We emphasize an important feature of the spectrum: it has more power at long wave- 
lengths because v — 3/2 > 0. This is recognized to be a consequence of the spinodal insta- 
bilities that result in the growth of long wavelength modes and therefore in more power for 
these modes. This seems to be a robust prediction of new inflationary scenarios in which the 
potential has negative second derivative in the region of field space that produces inflation. 

It is at this stage that we recognize the consistency of our approach for separating the 
composite effective zero mode from the small fluctuations. We have argued above that many 
more than 60 e-folds are required for consistency, and that the small fluctuations correspond 
to those modes that cross the horizon during the last 60 e-folds of the inflationary stage. 
For these modes H /k = e~ H o t *( k ) where t*(k) is the time since the beginning of inflation of 
horizon crossing of the mode with wavevector k. The scale that corresponds to the Hubble 
radius today Ao = 2Tr/k is the first to cross during the last 60 or so e-folds before the 
end of inflation. Smaller scales today will correspond to k > ko at the onset of inflation 
since they will cross the first horizon later and therefore will reenter earlier. The bound 
on \5k \ oc AT/T < 10~ 5 on these scales provides a lower bound on the number of e-folds 
required for these type of models to be consistent: 



N e > 60 + 



12 



ln(z/ 



v — 



(8.30) 



where we have written the total number of e-folds as N e = Hq t*(ko) + 60. This in turn can 
be translated into a bound on the coupling constant using the estimate given by eq. (|8.3|) . 

The four year COBE DMR Sky Map [^(J gives n w 1.2 ± 0.3 thus providing an upper 
bound on v 



< v- - < 0.05 
2 ~ 



(8.31) 



corresponding to h > 2.6. We then find that these values of ho and A ~ 10 12 — 10 14 
provide sufficient e-folds to satisfy the constraint for scalar density perturbations. 
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F. Tensor Metric Perturbations: 



The scalar field does not couple to the tensor (gravitational wave) modes directly, and 
the tensor perturbations are gauge invariant from the beginning. Their dynamical evolution 
is completely determined by the dynamics of the scale factor [37,41]. Having established 



numerically that the inflationary epoch is characterized by H / Hq <C 1 and that scales 
of cosmological interest cross the horizon during the stage in which this approximation is 
excellent, we can just borrow the known result for the power spectrum of gravitational waves 
produced during inflation extrapolated to the matter era p7| , |4l| 



Mk) « . (8.32) 

Thus the spectrum to this order is scale invariant (Harrison-Zeldovich) with an amplitude 
of the order m 4 /XMp l . Then, for values of m ss 10 12 — 10 14 Gev and A ss 10~ 12 — 10~ 14 
one finds that the amplitude is < 10~ 10 which is much smaller than the amplitude of scalar 
density perturbations. As usual the amplification of scalar perturbations is a consequence 
of the equation of state during the inflationary epoch. 



IX. CONCLUSIONS 

Since there are a number of articles in the literature treating related problems, it is 
useful to review the unique features of the present work. First, we have treated the problem 
dynamically, without using the effective potential (an equilibrium construct) to determine 
the evolution. Second, we have provided consistent non-perturbative calculations of the 
non-linear quantum field evolution to bring out some of the most relevant aspects of the late 
time behavior. In particular, we found that the quantum backreaction naturally inhibits 
catastrophic growth of fluctuations and provides a smooth transition to the late time regime 
in which the quantum fluctuations decay as the zero mode approaches its asymptotic state. 
Third, the dynamics studied obeys the constraint of covariant conservation of the energy 
momentum tensor. 

It can be argued that the inflationary paradigm as currently understood is one of the 
greatest applications of quantum field theory. The imprint of quantum mechanics is every- 
where, from the dynamics of the inflaton, to the generation of metric perturbations, through 
to the reheating of the universe. It is clear then that we need to understand the quantum 
mechanics of inflation in as deep a manner as possible so as to be able to understand what 
we are actually testing via the CMBR temperature anisotropies, say. 

What we have found in our work is that the quantum mechanics of inflation is extremely 
subtle. We now understand that it involves both non-equilibrium as well as non-perturbative 
dynamics and that what you start from may not be what you wind up with at the end! 

In particular, we see now that the correct interpretation of the non-perturbative growth 
of quantum fluctuations via spinodal decomposition is that the background zero mode must 
be redefined through the process of zero mode reassembly that we have discovered. When 
this is done (and only when!) we can interpret inflation in terms of the usual slow- roll 
approach with the now small quantum fluctuations around the redefined zero mode driving 
the generation of metric perturbations. 
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We have studied the non-equilibrium dynamics of a 'new inflation' scenario in a self- 
consistent, non-perturbative framework based on a large N expansion, including the dy- 
namics of the scale factor and backreaction of quantum fluctuations. Quantum fluctuations 
associated with superhorizon modes grow exponentially as a result of the spinodal insta- 
bilities and contribute to the energy momentum tensor in such a way as to end inflation 
consistently. 

Analytical and numerical estimates have been provided that establish the regime of 
validity of the classical approach. We find that these superhorizon modes re-assemble into 
an effective zero mode and unambiguously identify the composite field that can be used 
as an effective expectation value of the inflaton field whose classical dynamics drives the 
evolution of the scale factor. This identification also provides the initial condition for this 
effective zero mode. 

A consistent criterion is provided to extract small fluctuations that will contribute to 
cosmological perturbations from large non-perturbative spinodal fluctuations. This is an 
important ingredient for a consistent calculation and interpretation of cosmological pertur- 
bations. This criterion requires that the model must provide many more than 60 e-folds 
to identify the 'small perturbations' that give rise to scalar metric (curvature) perturba- 
tions. We then use this criterion combined with the gauge invariant approach to obtain the 
dynamics of the Bardeen variable and the spectrum for scalar perturbations. 

We find that during the inflationary epoch, superhorizon modes of the Bardeen potential 
grow exponentially in time reflecting the spinodal instabilities. These long-wavelength insta- 
bilities are manifest in the spectrum of scalar density perturbations and result in an index 
that is less than one, i.e. a 'red' power spectrum, providing more power at long wavelength. 
We argue that this red spectrum is a robust feature of potentials that lead to spinodal in- 
stabilities in the region in field space associated with inflation and can be interpreted as an 
imprint of the phase transition on the cosmological background. Tensor perturbations on 
the other hand, are not modified by these features, they have much smaller amplitude and 
a Harrison-Zeldovich spectrum. 

X. ACKNOWLEDGEMENTS: 

We thank F. Cao, D. Cormier, R. Holman, S. P. Kumar, J. Salgado, A. Singh, M. 
Srednicki all of whom collaborated at different stages on the works reviewed here. We thank 
J. Baacke, C. Destri, A. Dolgov, E. Kolb, E. Weinberg for conversations and discussions. 
D. B. thanks the N.S.F for partial support through the grant awards: PHY-9605186 and 
INT-98 15064 and LPTHE for warm hospitality. We thank the CNRS-NSF cooperation 
programme for partial support. 



57 



REFERENCES 



Laboratoire Associe au CNRS, UMR7589. 
[1] A. H. Guth, Phys. Rev. D23, 347 (1981). 

[2] For thorough reviews of standard and inflationary cosmology see: E. W. Kolb and M. 
S. Turner, The Early Universe (Addison Wesley, Redwood City, C.A. 1990). A. Linde, 
Particle Physics and Inflationary Cosmology, (Harwood 1990) and ref. 0. 

[3] For more recent reviews see: M. S. Turner, astro-ph-9703197; astro-ph-9703196; astro- 
ph-9703174; astro-ph-9703161; astro-ph-9704062; astro-ph-9704024. A. Linde, in Cur- 
rent Topics in Astrofundamental Physics, Proceedings of the Chalonge Erice School, N. 
Sanchez and A. Zichichi Editors, Nato ASI series C, vol. 467, 1995, Kluwer Acad. Publ. 
A. R. Liddle, astro-ph-96 12093, Lectures at the Casablanca School Morocco, 1996. 

[4] G. Smoot, in the Proceedings of the Vth. Erice Chalonge School on Astrofundamental 
Physics, p. 407-484, N. Sanchez and A. Zichichi eds., World Scientific, 1997. 
A. R. Liddle and D. H. Lyth, Phys. Rep. 231, 1 (1993). 

A. E. Lange et al. , |astro-ph/000500l , S. Hanany et al., |astro-ph/0005T2l , P. de 

Bernardis et al., Nature 404 (2000) 955. 
[5] for reviews of inflation, see R. Brandenberger, Rev. of Mod. Phys. 57, 1 (1985); Int. J. 

Mod. Phys. A2, 77 (1987) and ref. @. 
[6] D. Boyanovsky, D. Cormier, H. J. de Vega, R. Holman, S. P. Kumar, 'Out of equilibrium 

fields in inflationary dynamics, density fluctuations', in Proceedings of the NATO ASI 

'Current Topics in Astrofundamental Physics: Primordial Cosmology', edited by N. 

Sanchez and A. Zichichi, vol C511, Kluwer Dordrecht, 1998. 
[7] D. Boyanovsky and H. J. de Vega, Phys. Rev. D47, 2343 (1993). 
[8] D. Boyanovsky, H. J. de Vega, R. Holman, D.-S. Lee and A. Singh, 

Phys. Rev. D51, 4419 (1995); 

D. Boyanovsky, M. D'Attanasio, H. J. de Vega, R. Holman and D. S. Lee, 
Phys. Rev. D52, 6805 (1995); 

For reviews see, D. Boyanovsky, H. J. de Vega and R. Holman, in the Proceedings of 
the Second Paris Cosmology Colloquium, Observatoire de Paris, June 1994, p. 127- 
215, H. J. de Vega and N. Sanchez Editors, World Scientific, 1995; D. Boyanovsky, M. 
D'Attanasio, H. J. de Vega, R. Holman and D.-S. Lee, 'New aspects of reheating', in the 
Proceedings of the Erice Chalonge School, 'String Gravity and Physics at the Planck 
Energy Scale', NATO ASI, N. Sanchez and A. Zichichi Editors, Kluwer 1996, p. 451-492. 
[9] D. Boyanovsky, H.J. de Vega, R. Holman, J.F.J. Salgado, Phys. Rev. D54, 7570 (1996). 
D. Boyanovsky, H. J. de Vega and R. Holman in the Proceedings of the Vth. Erice 
Chalonge School on Astrofundamental Physics, p. 183-270, N. Sanchez and A. Zichichi 
eds., World Scientific, 1997. 

D. Boyanovsky, C. Destri, H.J. de Vega, R. Holman and J.F.J. Salgado, 

Phys. Rev. D57, 7388 (1998). 
[10] D. Boyanovsky, H. J. de Vega, R. Holman and J. Salgado, 

Phys. Rev. D59, 125009 (1999). 
[11] D. Boyanovsky, H. J. de Vega, and R. Holman, Phys. Rev. D49, 2769 (1994). 
[12] D. Boyanovsky, D. Cormier, H. J. de Vega, R. Holman, A. Singh, M. Srednicki, 

Phys. Rev. D56, 1939 (1997). 
[13] D. Boyanovsky, D. Cormier, H. J. de Vega and R. Holman, 



58 



Phys. Rev. D55, 3373 (1997). 
[14] D. Boyanovsky, D. Cormier, H. J. de Vega, R. Holman and S. P. Kumar, Phys. Rev. 

D57, 2166 (1998). 
[15] D. Boyanovsky, M. D'Attanasio, H. J. de Vega and R. Holman, 

Phys. Rev. D54, 1748 (1996), and references therein. 
[16] D. Boyanovsky and H. J. de Vega, Phys. Rev. D61, 105014 (2000). 
[17] D.T. Son, Phys. Rev. D54, 3745 (1996); |hep-ph/960l3T7| . 

[18] D.I. Kaiser, Phys. Rev D53, 1776 (1996), D56, 706 (1997) and D57, 702 (1998). 

[19] J. Schwinger, J. Math. Phys. 2, 407 (1961); P. M. Bakshi and K. T. Mahanthappa, J. 

Math. Phys. 4, 1 (1963); ibid, 12; L. V. Keldysh, Sov. Phys. JETP 20, 1018 (1965); A. 

Niemi and G. Semenoff, Ann. of Phys. (N.Y.) 152, 105 (1984); Nucl. Phys. B [FS10], 

181 (1984); E. Calzetta, Ann. of Phys. (N.Y.) 190, 32 (1989); R. D. Jordan, Phys. Rev. 

D33, 444 (1986); N. P. Landsman and C. G. van Weert, Phys. Rep. 145, 141 (1987); 

R. L. Kobes and K. L. Kowalski, Phys. Rev. D34, 513 (1986); R. L. Kobes, G. W. 

Semenoff and N. Weiss, Z. Phys. C29, 371 (1985). 
[20] For non-equilibrium methods in cosmology see, for example: E. Calzetta and B.-L. Hu, 

Phys. Rev. D35, 495 (1988); ibid D37, 2838 (1988); J. P. Paz, Phys. Rev. D41, 1054 

(1990); ibid D42, 529 (1990); B.-L. Hu in Bannf/Cap Workshop on TFT, ed. by F.C. 

Khanna, R. Kobes, G. Kunstatter, H. Umezawa, World Scientific, Singapore, 1994, 

p. 309 and in the Proceedings of the Second Paris Cosmology Colloquium, Observatoire 

de Paris, edited by H. J. de Vega and N. Sanchez (World Scientific, Singapore, 1995), 

p. Ill and references therein. 
[21] H. Leutwyler and S. Mallik, Ann. of Phys. (N.Y.) 205, 1 (1991). 
[22] O. Eboli, R. Jackiw and S.-Y. Pi, Phys. Rev. D37, 3557 (1988); 

M.Samiullah, O. Eboli and S-Y. Pi, Phys. Rev. D44, 2335 (1991). 

J. Guven, B. Liebermann and C. Hill, Phys. Rev. D39, 438 (1989). 
[23] N. D. Birrell and P.C.W. Davies, Quantum fields in curved space (Cambridge Univ. 

Press, Cambridge, 1986). 
[24] J. Baacke, K. Heitmann and C. Patzold, Phys. Rev. D55, 2320 (1997), D56, 6556 

(1997), D57, 6398 and 6406(1998). 
[25] M. Abramowitz and I.E. Stegun (eds.), Handbook of Mathematical Functions (National 

Bureau of Standards, Washington, D.C., 1972), chapter 13. 
[26] N. Turok, Ue-Li Pen, U. Seljak, Phys. Rev. D58, 023506 (1998). A. Albrecht, R. A. 

Battye, J. Robinson, Phys. Rev. Lett. 79, 4736 (1997). B. Allen, R. R. Caldwell, S. 

Dodelson, L. Knox, E. P. S. Shellard, A. Stebbins, Phys. Rev. Lett. 79, 2624 (1997). R. 

Durrer, |astro-ph/ 0003364 
[27] D. H. Lyth, hep-ph-9609431 (1996). S. Dodelson, W. H. Kinney and E. W. Kolb, Phys. 

Rev. D56, 3207 (1997). 
[28] A.D. Linde, Phys. Lett. B116, 335 (1982). A. Vilenkin and L. H. Ford, Phys. Rev. D26, 

1231 (1982). A. Vilenkin, Nucl. Phys. B226, 504 (1983); Nucl. Phys. B226, 527 (1986). 
[29] A. Vilenkin, Phys. Lett. B115, 91 (1982). 

[30] P. J. Steinhardt and M. S. Turner, Phys. Rev. D29, 2162, (1984). 
[31] A. Guth and S-Y. Pi, Phys. Rev. D32, 1899 (1985). 

[32] For non-equilibrium methods in different contexts see for example: F. Cooper, J. M . 
Eisenberg, Y. Kluger, E. Mottola, B. Svetitsky, Phys. Rev. Lett. 67, 2427 (1991); F. 



59 



Cooper, J. M. Eisenberg, Y, Kluger, E. Mottola, B. Svetitsky, Phys. Rev. D48, 190 
(1993). 

[33] F. Cooper and E. Mottola, Mod. Phys. Lett. A 2, 635 (1987); F. Cooper, S. Habib, Y. 
Kluger, E. Mottola, J. P. Paz, P. R. Anderson, Phys. Rev. D50, 2848 (1994). F. Cooper, 
S.-Y. Pi and P. N. Standoff, Phys. Rev. D34, 3831 (1986). F. Cooper and E. Mottola, 
Phys. Rev. D36, 3114 (1987). 

F. Cooper, Y. Kluger, E. Mottola, J. P. Paz, Phys. Rev. D51, 2377 (1995). 
[34] S. A. Ramsey, B. L. Hu, Phys. Rev. D56, 678 (1997). 

[35] D. Boyanovsky, H. J. de Vega and R. Holman, Phys. Rev. D 51, 734 (1995). 
D. Boyanovsky, H. J. de Vega, R. Holman and S. Prem Kumar, 
Phys. Rev. D56, 3929 and 5233 (1997). 
The last reference under [[33]. 

J. Lesgourgues, D. Polarski and A. A. Starobinsky, Nucl. Phys. B497, 479 (1997). 
[36] L. P. Grishchuk, Phys. Rev. D 45, 4717 (1992). 

[37] V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rep. 215, 293 (1992). 
[38] J. Bardeen, Phys. Rev. D22, 1882 (1980). 

[39] R. R. Caldwell, Class. Quant. Grav. 13, 2437 (1995). J. Martin and D. J. Schwarz, Phys. 
Rev. D57 (1998) 3302. 

[40] K. M. Gorski, A. J. Banday, C. L. Bennett, G. Hinshaw, A. Kogut and G. F. Smoot, 
astro-ph-9601063. 

[41] L. P. Grishchuk, Phys. Rev. D. 52, 5549, (1995); Proceedings of the Erice Chalonge 
School NATO ASI on 'String Gravity and Physics at the Planck Scale', Ed. N. Sanchez 
and A. Zichichi, (Kluwer, 1996), p. 369; Phys. Rev. D53 (1996) 6784; Proceedings of 
the NATO ASI on 'Current Topics in Astrofundamental Physics', Ed. N. Sanchez and 
A. Zichichi (Kluwer, 1995), p. 205. 
Rev. Mod. Phys. 69, 373, (1997). 

[42] I. Zlatev, G. Huey and P. J. Steinhardt, Phys. Rev. D 57, 2152 (1998). 

[43] F. Cooper, S. Habib, Y. Kluger, E. Mottola, J. P. Paz, and P. R. Anderson, Phys. Rev. 
D50, 2848 (1994). 

[44] J. Avan and H. J. de Vega, Nucl. Phys. B224, 61 (1983). 
M. B. Halpern, Nucl. Phys. B 173, 504 (1980). 

[45] D. Boyanovsky, H. J. de Vega, R. Holman, S. Prem Kumar and R. D. Pisarski, Phys. 
Rev. D57, 3653 (1998). 

F. J. Cao and H. J. de Vega, |hep-ph/99ll537. 



60 



